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Abstract 

A systematic derivation is presented of the twist-2 anomalous dimensions of the general 
quark and gluon light-ray operators in the generalized Bjorken region in leading order 
both for unpolarized and polarized scattering. Various representations of the anomalous 
dimensions are derived in the non-local and local light cone expansion and their properties 
are discussed in detail. Evolution equations for these operators are derived using different 
representations. General two- and single-variable evolution equations are presented for the 
expectation values of these operators for non-forward scattering. The Compton amplitude 
is calculated in terms of these distribution amplitudes. In the limit of forward scattering a 
new derivation of the integral relations between the twist-2 contributions to the structure 
functions is given. Special limiting cases which are derived from the general relations are 
discussed, as the forward case, near-forward scattering, and vacuum-meson transition. 
Solutions of the two-variable evolution equations for non-forward scattering are presented. 



1 Introduction 



The Compton amplitude for the scattering of a virtual photon off a hadron 

7* + Pi^l'* +P2 (1.1) 

provides one of the basic tools to understand the short-distance behavior of the nucleon and to 
test Quantum Chromodynamics (QCD) at large space-like virtualities. The Compton amplitude 
for the general case of non-forward scattering is given by 



Here. 



T^(p+,p.,q)=i / (fxe tqx {p 2 ,S 2 \T(J fl (x/2)J v (-x/2))\p 1 ,S 1 ) . (1.2) 



p+ = P2+P1, V- =Vi -Pi = 9i - 92, (1-3) 
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0(91 + 92), Pi + qi=P2 + q2, (1.4) 



where q\ (q 2 ) and p\ (p 2 ) denote the four-momenta of the incoming (outgoing) photon and 
hadron, respectively, and Si, 5*2 are the spins of the initial- and final-state hadron. The gener- 
alized Bjorken region is defined by the conditions 

v = qp + — > 00, Q 2 = —q 2 — > 00 , (1.5) 

keeping the variables 

£ = > »7 = = — ^ (I- 6 ) 

9P+ 9P+ 2z/ 

fixed. Relations between these kinematic variables are given in Appendix |A] For = one 
obtains rj = — £. In distinction to deep inelastic forward scattering, the scaling variable £ is 
not restricted to the interval < £ < 1. For forward scattering q = q± = q 2 ,p+ = 2p,r] = 0, 
£ - ► ^Bj = Q 2 / (2p9), holds. In principle the results presented below can be analytically continued 
to time-like processes with two space-like initial state photons as a necessary condition. 

The (renormalized) time-ordered product in Eq. (|1.2j ) can be represented in terms of the 
operator product expansion. Unlike the local operator product expansion [I], which is widely 
used in the case of forward scattering, the non-local operator product expansion [2-5] , to which 
we refer in the present paper, leads to compact expressions for coefficient functions and oper- 
ators in the non-forward case. As was shown in Ref. the expressions of the local operator 
product expansions can be obtained by a Taylor expansion of the non-local operator product 
expansion (cf. also 0, §). In lowest order in the coupling constant the expansion contains only 
quark operators with two external legs. These operators can be decomposed into operators of 
different twist. In the present paper we will consider the contribution of the twist-2 operators 
only and calculate the anomalous dimensions which emerge in the corresponding renormalization 
group equations. 

Our goal is to demonstrate that quite different processes, such as deep inelastic forward scat- 
tering, the different channels in deeply virtual Compton scattering, i.e. non-forward scattering 
at a general kinematics, cf. [7-9] for early studies, and vacuum-meson transitions contain as an 
essential input the same light-ray operators of twist-2. Therefore the scaling violations of these 
processes are described by the same anomalous dimensions. The non-perturbative partition 
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functions are obtained as the Fourier transforms of the corresponding matrix elements of these 
operators. With the help of the Fourier representations we derive a new representation for the 
Compton scattering amplitude in the generalized Bjorken region which extends the representa- 
tion known for forward scattering. Moreover, the partition functions as the non-perturbative 
input distributions are related to the Fourier transforms of their matrix elements. Therefore, the 
evolution equations result directly from the renormalization group equation of the light-ray op- 
erators involved. Consequently, all the evolution kernels can be determined from the anomalous 
dimensions of the corresponding operators. As special cases one obtains the evolution equations 
for forward scattering and for the case of vacuum-meson transitions. Furthermore the equiv- 
alence to other representations in the recent literature is shown, in which different kinematic 
variables have been used. 

In the generalized Bjorken region the distribution amplitudes for the twist-2 contributions 
depend on two scaling variables unlike the case of forward scattering. One of these parameters 
is related to the sum, the other to the difference of the initial and final state nucleon momenta. 
In [|TI], [TT| a solution of the evolution equation for the non-singlet case was obtained. In the 



literature mostly the one-parametric equations have been discussed so far by imposing further 



conditions on the momenta. The solution for the singlet case was given in Ref. |12| for the 
two-variable equations. 

The paper is organized as follows. In Section 2 we derive the non-local light cone expansion 
of the Compton amplitude for non-forward scattering in lowest order. The twist decomposition 
of the non-local operators is performed in Section 3. In Section 4 the non-forward matrix 
elements of the non-local operators are calculated and related to the corresponding distribution 
amplitudes. The double and single-variable distribution amplitudes contributing to the Compton 
amplitude are discussed in Section 5, where also a derivation of special relations holding in the 
case of forward scattering off unpolarized and polarized nucleons is presented. In Section 6 a 
detailed derivation of the anomalous dimensions of the quark and gluon operators contributing 
at the leading twist level to the non-forward Compton amplitude are given in 0(a s ). The 
evolution kernels are presented in different representations both for unpolarized and polarized 
scattering. The corresponding evolution equations are discussed in Section 7, both for the case of 
two and one-variable evolution equations. Here also the kernels for the single-variable equations 
are given. Special cases of evolution equations into which the general equations transform in a 
series of kinematic limits are discussed in Section 8, including the cases of forward scattering 
and vacuum-meson transitions. The solutions of the general two-variable evolution equations 
are provided in Section 9. Section 10 contains the conclusions. In the Appendices details of the 
calculation are presented. A brief summary of part of our results has already appeared In 
this paper the derivation is presented in detail. 



2 The non- local light— cone expansion in Born approxi- 
mation 

The Compton amplitude will be evaluated using the non-local light cone expansion (LCE) in 
terms of bi-local light-ray operators. This expansion is a summed-up version of the local LCE. 
The latter can be obtained performing a Taylor expansion of the former representation. The 
rigorous proof of the non-local LCE is rather complicated. Here, following the prescription of 
Refs. ||, H , it will be introduced heuristically with some comments on the essential steps. 
Firstly, one considers the renormalized (R) time ordered (T) product of two electromagnetic 
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currents 

RT{J,{x/2)J v {-x/2)S), (2.1) 

which in perturbative QCD is represented on the Hilbert space of free fields. The electromagnetic 
current reads 

J,(x) = ^(x) 7At A en V(x) , 
where ip(x) are the quark fields and 

Aem = K A ? + 73 A ?) for SU ^ 



flavor 



S denotes the renormalized S-matrix. 

In lowest order in the coupling constant we have S = 1 and, expressing the T-product in 
terms of normal products, we obtain 

RT{J,{x/2)J v {-x/2)S) = t:^(x/2){S, upa -te^ pal5 }r(*d p x D%x))c a \ a f ^(-x/2) : 

- \(x/2, fi) <- (-x/2, u)] (2.2) 
+ higher order terms , 

where 

Spupcr — S^pi/u QixvQpa Qp,pQvcr Qp/rQpv 

and Sp Upa denotes the Levi-Civita symbol. The S77(3)fl av0 r vector 

with e the electric charge, determines the flavor content and D c (x) = (— i/4:7r 2 )(x 2 — ie)~ l is the 
free scalar propagator. Because this expression is restricted to leading order no phase factors 
occur. 

Secondly, since the non-local LCE of composite operators appears as a (multiple) integral 
representation we hint at it by introducing two auxiliary /t-integrations [||] [] 

I C a X p r+oo r+oo 

RT(Jp(x/2)J u (-x/2)S) = — — — / d«i / dK 2 (2.3) 

Z7T yX — IB) J—oo J—oo 
X {§ ^{Kix)X a f ^ll){K2x) - ^(K 2 x)X a f J a l/j(K 1 x)\ (-S^pa) A- («i , K 2 ) 

+~ i)(K l x)^\ a p a i)(K 2 x) + i)(K 2 x)^\ a p a i)(K l x) ie^A+^K^j 

+ ... , 



where 



A ± (ki, k 2 ) = 5(aci - \)5(k 2 + \)± 5(k 2 - + h 



The next step consists in the projection of the above expressions onto the light-cone by 
replacing x — > x with x 2 = 0, leaving x 2 unchanged. The light-like vector 



x = x + — 

c 



vW) 2 - x 2 c 2 - «)) 



1 Here, and in the following, the normal product symbols for the relevant operator products will be omitted. 
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is related to x and a fixed non-null subsidiary four-vector ( the dependence of which drops out in 
leading order expressions. It has to be noted that because the light-cone expansion is considered 
for the operators and not, as often used, for their matrix elements, a light-like reference vector 
has necessarily to be introduced in configuration space.Q 

In general, also non-leading contributions related to higher order terms of the S'-matrix and, 
in certain gauges, also of the phases factors U(x/2, —x/2) have to be taken into account. Then 
the coefficient functions like (c a /27r 2 )(x 2 — ie) 2 A±(K,i, K2) depend non-trivially on /tj. In this 
way one arrives at a pre-form of the non-local light-cone expansion being given by a sum over 
the singular coefficient functions C a (x 2 , Hi) times the light-ray operators O a (Kix), 

J2C a {x 2 ,Ki) X O a {KiX) . 
a 

It should be remarked that even if the perturbative expansion of the renormalized time ordered 
product of the two currents has been correctly performed, there is no reason to assume that the 
appearing light-ray operators are the renormalized ones: The expansions resulting from Eq. ( |2.3|) 
are not yet the true light-cone expansions. 

A complete proof which shows that the free field operators appearing in Eq. (|2.3| ) have to be 
substituted by renormalized light-ray operators is much more involved ^ |3| . Firstly, an operator 
identity for the product of two composite operators, e.g., the two currents, will be proven which 
holds on the whole Hilbert space. Then, introducing a new light-cone adapted renormalization 
procedure R it can be shown that the perturbative functional of the renormalized product of 
two currents may be split into an asymptotically relevant part - the light-cone expansion up to 
a definite order in the light-cone singularity - and a well-defined remainder being less singular: 

RT (J^(x/2)J u (-x/2)S)^ J + ^d 2 KC T {x 2 ,^^ 2 )RT{p T {K l i,K 2 x)S) + higher order terms, (2.4) 
where the free light -ray operators with a specified T-structure are given by 

O r (^, n 2 x) = / ( S 4 ( ^)4 e™ +J ™ : V^i)rV(p 2 ) : ■ (2.5) 

Here, the coefficient functions and the corresponding light-ray operators appear quite naturally 
as renormalized ones. In the following, for brevity, we usually do not explicitly indicate that 
they are renormalized quantities. 

Let us point to the fact that, in principle, within this formalism the non-local as well as the 
usual local light-cone expansion may be obtained (see also |J). The only difference consists in 
the following: In the non-local LCE, as an intermediate step to obtain the final representation 
Eq. Q2.4| ), a Fourier transformation of the renormalized, perturbatively determined coefficient 
functions Fr(x 2 ,xpi, /1 2 ) which multiply the various normal products of free operators in mo- 
mentum space, like : ip(pi)Tip(p2) : in Eq. ( |2.5| ), has been performed. Primarily, these functions 
are given as 

F r (x 2 ,x Pi ,fi 2 ) = J + ' d 2 Ke tK1 * pl+iK2 * P2 C r (x 2 ,K;v 2 ). (2.6) 

2 x and its pendant x* with xx* — 1 correspond to the light -like vectors n and p of dimension mass -1 and mass, 
respectively, usually introduced in deep inelastic scattering according to = + (M 2 /2)n fl [Q. Furthermore, 
x and x* define in a sense a 'comoving' reference system on the light cone. 
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In addition, by the ct-representation it can be proven that they are entire analytic functions 
with respect to the variables xpi M [5] and therefore the range of «j is restricted to 

- 1 < «i < +1 . (2.7) 

On the other hand, in the local LCE one applies a Taylor expansion of F-p with respect to xpi 
which leads to a (twofold) infinite series over the local operators. 

These local and non-local operators are connected by the following relations: 



T (kiX, K 2 x) 



With the convention 



{xBy 



/ y „ i „ I ninjv /i 



Kl=K2=0 



(2.8) 
(2.9) 



and choosing the axial gauge, x^A^x) = 0, we obtain 

o r nin2 (x) = m(d £) n T(s orm . (2.10) 

Summarizing the foregoing procedure and introducing, for later convenience, the variables 

k± — 2(^2 i «i) with = «+ (2-11) 

the general expression will be, taking into account the explicit form of S^ upa , 

RT(J^x/2)J u (-x/2)S) ~IJdk (2.12) 

+ iC aj5 (x 2 , k + , k_, /i 2 )e M ^ p<T a;p05 )(7 (K + x, /i 2 ) + higher order terms, 
where the measure 



£>K = dKidK 2 0(l-Ki)0(l + Ki)0(l -K 2 )0(1 + K 2 ) 

= 2ck + ck_0(l + k + + k_)0(1 + k+ - «_)0(1 - k + + «_)0(1 - k. 



(2-13) 



has been introduced. Here C a rg) are the renormalized coefficient functions which in Born ap- 
proximation read 



C a {x , k+, kJ) 

C a ,5( x j K +j K -) 



1 



27r2 (x2 _ i£ /(* + ) [*(*- - 1) - *(«- + 1) 

<J(«+) [<?(«_ - I) +5(k_ + §) 



27r 2 (x 2 — is) 2 



and 0? 5 \ CT are the renormalized (anti)symmetric bi-local LC-operators (] 



'(5) 



ip(Kix)\ a p a U (ki, k 2 )^/j(k 2 x) - vjj(K 2 x)\ a f 'y (T U (« 2 , «i)'0(«l5) 



5. 



(2.14) 
(2.15) 

(2.16) 



0^ a (Ki,K 2 ) = \RT V(Kix)A^757 CT f/(Ki,/t 2 )^(K 2 x) + ip(K 2 x)X° f 'j 5 'y (7 U(K 2 ,K 1 )i;(K 1 x) S, (2.17) 



3 In the following (k\, K2) and (k+, K-) are chosen interchangeably, i.e. 0(k + x, K—x) is written instead of the 
more lengthly expression 0((k+ — K-)x, (k + + K—)x). Furthermore, very often x and the renormalization scale 
[i 2 will be suppressed in the arguments of the operators and related quantities for brevity. 
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respectively, with the phase factors 



U(ki, k 2 ) = U(kix, k 2 x) — Vexjp |— ig J dr x M v4 M (ra;)| 



(2.18) 



where V denotes the path ordering, g is the strong coupling constant and = A a ^t a is the gluon 
field with t a being the generators of S77(3) co i or in the fundamental representation spanned by 
the quark fields. In the explicit expressions we have indicated by RT[0]S that these operators 
are renormalized. 

The general expression Eq. (|2.12|) contains a flavor singlet as well as the flavor non-singlet 
(NS) parts. In the case of flavor non-singlet operators, i.e., for \°j 7^ 1, the consideration of the 
quark operators Eqs. fl2.16| ) and (|2.17| ) is sufficient. However, the singlet operators (A° = 1), in 
the following denoted by O q , mix with certain gluon operators which will be denoted by G . 
In addition, if we would consider the Compton amplitude beyond leading order also the gluon 
operators would appear in the LCE. 

The gluon operators related to the quark operators Eqs. (|2.16|) and ( p. 17] ) read 



where, 



FZ x ( Kl x)u ab ( Kl , K 2 )F b vX {K 2 ~x) + f; a 



„ ■■ , -„ K K 2 x)U ah {K 2 ,K X )F h vX (K 2 X 

\RT [Fffaxprtfa, k 2 )F* x (k 2 x) - F,? A 



ah J 
K 2 x)U ab < 



K 2 , 



S, (2.19) 
S, (2.20) 



p — p a + a 

1 )1U — 1 fj, u L 



and F, 



-£ F pa 



is the gluon field strength and the dual field strength, respectively. U ab (K 1 ,K 2 ) is the phase 
factor ( 2.18|) in the adjoint representation. To obtain exact correspondence to the above quark 
operators the gluon operators have to be multiplied by ^(g a . tJ 'x u + gjx^ 1 ). 



3 Twist decomposition of non— local operators 

The non-local operators Eqs. ( |2.16 , |2.17| , [2.19| , |2.20|) contain contributions of different twist. 



Here we use the notion of twist in its original definition |T5[ as geometric twist, i.e. canonical 
dimension minus spin of the corresponding operator. Let us now extract the leading twist 
parts, which are relevant for the discussion of virtual Compton scattering. This concerns the 
contributions of twist two and three, the latter being especially relevant for the antisymmetric 
part of the scattering amplitude which concerns the polarized scattering. 

We consider first the quark operators O a a and 0^ a , which contain contributions of twist-2, 3 
and 4: 

0a £^a,twist2 _|_ twist3 _|_ ^a,twist4 (3 1) 

Oa ^~ja,twist2 , ^a, twist3 , / ^a,twist4 /q r\\ 

5,cr — U 5,cr + U 5,cr + U 5,cr ■ 

In principle also higher twist contributions appear, but they vanish on the light-cone since they 
are proportional to nonvanishing integer powers of x 2 . 

To obtain the bi-local operators (|3.1|) of definite twist we consider the centered operators 



ip(—Kx)\J'y a U(—Kx,Kx)ilj(Kx) (3.3) 
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with arbitrary values of x 7^ x, i.e., we consider the operators as not restricted to the light-cone. 
They are represented through the Taylor expansion Eq. ( |2.8| ) by a series of local tensor operators 
of order n = ri\ + n 2 which, in the axial gauge, may be written as 



Here we used the abbreviation d=d — d- By purely group theoretical methods these local 
operators are decomposed into traceless tensor operators of definite dimension and spin carrying 
an irreducible representation of the Lorentz group. The symmetry behavior of these tensors is 
characterized by a few special classes of Young tableaux'. As is well known, the twist-2 terms 
correspond to the total symmetric representation. Finally, the local operators of definite twist 
are resummed to give non-local harmonic tensor operators of definite twist. 

Technically, these non-local harmonic tensor operators are obtained from the initial ones 
by first projecting onto traceless operators and then applying appropriate differential operators 
related to the corresponding symmetry class. Finally, the resulting expressions are projected 
onto the light-ray x. In principle this method can also be applied for arbitrary values of ftj. 

This approach has been worked out in Refs. \TE, |T7[ for various light-ray operators. Here, we 
refer only to those representations which are relevant for the subsequent considerations. First 
the projections onto the traceless operators are defined. For the non-local scalar and vector 
quark operators Eqs. ( |2.16|) and ( [2.171) , respectively, they are given by 



O 



q, traceless / 



^ \^p{— K,x)(x / y)ip(K,x) — ip(K,x)(x r y)ip(—Kx) 



(3.4) 



+ £ / dt 



k=l 



1 -t 



k-1 



-X 



k\(k- 1)! 



X f [^( _ Ktx)(x r y)lj)(Ktx) — 1p(Ktx)(x^)'lp(— Ktx) 



Qq, traceless ^ = | ^-^^{kx) ~ ^(kx)^(-Kx) 

CO 



(3.5) 



+ U 7 



1 dt n - t\ k -W -x 2 \ n k 



+ 



t J \ 4 J k\(k-l)\ 

2 — Ktx)"f a 1p(Ktx) — tfj(K,tx)'~f a llj( — Ktx) 

2 [(^XpD + x 2 d a d p ) - 2x a d p (xd) 



x 



^ <■' - 1 '1 - t\ k ( -x 2 \ k □ " 



x ]T / dr I dt 



k=0 











k\k\ 



X § \^^~ KT tx) , y P ll ! {^Ttx) — 1p(K,Ttx)'~f P ll ! ( — K,Ttx) 

Here the phase factors and the generators \J of the flavor group have been suppressed and the 
operators generically are indicated by q. The sum over terms containing the d'Alembertian are 
the subtractions of the traces. The representation for O q ' traccloss ; Eq. ( |3.4|) , has been obtained 
already earlier in Ref. |T^|. Taking the limit x — > x in Eq. ( |3.4| ) the additional terms vanish 
because of x 2 = 0. Therefore the operator x a O^{ of twist two: 



o 



q, twist2 / 



o 



q, traceless 



(3.6) 



8 



Because of their tracelessness the above operators fulfill the following relations 

DO 9 ' traceless (-ra, kx) = 0, d a O q a > traceless (-ra, kx) = 0, 00% tracclcss (-ra, kx) = 0. (3.7) 

which proves that they are harmonic scalar and vector operators, respectively. The differential 
operators, to be applied to O 9, tracclcss and O q p ' traceless in order to generate the correct symmetry 
behavior, are d a for twist-2 and (g ap (xd)—x u d p ) for twist-3, respectively. Note that taking these 
derivatives before carrying out the light-cone limit leads to non-trivial contributions which result 
from the trace terms. 

The explicit computations straightforwardly lead to the following expressions for the twist-2, 
twist-3 and twist-4 light-ray vector operators |16|, |I7|: 



O 



g, twist 2 / 



Q q twist3 /_ ~ ^ 



o 



q, twist4 / 



dr d a O q 



dr 



traceless 

d a + i(lnr)x CT D 



9 ( — KTX, KTx) 



IX=X ' 



I dr (g ap xd - x^dp) O™ 

<j 

/ dr [{g ap xd - x a d t 

J 



traceless I KTX, KTX j 



\x=x ' 



- x a ((l + 21nr)<9 p + (lnr)xpD)] 9 ' p (-ktx, ktx) 



dr 



1 + \nT)d p + \{\UT)X P U 



o 



"PI 



-KTX, KTX) 



\x=x ' 



x=x ' 



(3.8) 



(3.9) 



(3.10) 



where the last expression is obtained using Eq. (|3.1| ). Note that the operators appearing at 
the r.h.s. of Eqs. ( |3.8| - |3.10| ) are to be taken with the phase factors before carrying out the 
derivatives. 

Important properties of these operators are 



x 



a ^g,twist2 O 9 ' 



twist 2 



~a Qq twistS = Q 



~™twist4 = 



(3.11) 



■\q, twist2 ^(j, twist3 



Analogous representations hold for 5 a , ^ 5(J 



On, twisw i y oc/,twist4 l • l l -i i - 

k'„ and Of a , which have similar properties. 



Note that the operators ( |3.8| - pTTU| ) include the trace terms which are proportional to x a . In fact 
they are again summed-up local operators fulfilling Eq. (|2.8|) . 

Let us now consider the gluon operators Eqs. ( |2.19| ) and Q2.2U| ). Again they contain a twist-2, 
twist-3 and twist-4 part, related to the traceless symmetric, the antisymmetric and the trace 
part of the tensor 0%^ , respectively. Since we are concerned with the scalar and vector gluon 
operators, 



X X U (5)fj,u 



and \{g^ + g^)0% )ian 



only the symmetric twist-2 part will be considered here. In the same manner as in the case of 
quark operators we first define the projection onto the centered traceless operator. It is given by 



O 



G, tracclcss / 



2 x^Xi, 



kx)F« (kx) + F* x (kx)F* : 



+ E 

k=l 



dtt 



—x 



2\* 



-KX 



(3.12) 



k\(k- 1)! 



"\/ /-yi /-yi 



F a p \-Ktx)F* x (Ktx) + F«\Ktx)F« x (-Ktx) 
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where again the phase factors have been suppressed. It fulfills the relation 



no 



G, traceless , 



o. 



(3.13) 



In comparison with Eq. ( |3.4| ) only an additional factor of t appears under the integral. This 
results from the fact that the operator (2.iy) relative to ( [2.16 ) contains one more factor of x not 
being accompanied by the parameter k. The scalar and vector twist-2 gluon light-ray operators 
are now given by 



G ' twist2 (-KX,Kx) 

0?> twiSt2 (-KX,Kx) 



qG, traceless ^ = G (~KX, Kx), 

f 1 drrda G ' traceless ( 

J 



drr 



G ( — KTX, KTX] 



(3.14) 



(3.15) 



Here we gave for simplicity the expressions for operators with the special arguments K\ = 
—k, k>2 = However, any of the Eqs. ( |3.4| - |3.15| ) remains true if general arguments are 
chosen. As an example the singlet vector operators of twist-2 read 



O 



A, twist2 / 



[K\X, K2X) 



drr 



d A -\ 



o- 



[K\TX, K2TX, 



(3.16) 



where A = (q, G) and d q — 1, do — 2. In addition, let us emphasize that in the various integrals 
defining operators of definite twist the variables have to be scaled but not the coordinates x 
of the corresponding (original) operators. 

The general form of harmonic operators of definite twist is necessary for an investigation of 
their matrix elements whereas the special centered form is sufficient to determine their anomalous 
dimensions, as will be shown in Sections 4 and 6, respectively. Inserting the result obtained for 
the twist-2 terms, Eq. (|3.8|), into Eq. fl2.12| ) we obtain an asymptotic representation of the 
Compton scattering amplitude Eq. (|1.2| ) which is valid in the generalized Bjorken region, as will 
be shown in Section 5 in detail. Concerning the renormalization of the twist-2 singlet operators 
and the evolution equations of the corresponding distribution amplitudes also the representations 
Eq. Q3.14|) and (|3.15|) necessarily come into play. 



4 Matrix elements of non-local operators 

The Compton scattering amplitude is obtained by taking the matrix elements of the operator, 
Eq. ( |2.12| ), between the hadron states \pi, Si). Thereby the matrix elements of the light-ray oper- 
ators Eqs. ( 2.16 ) and (|2.17|) appear as non-perturbative inputs. These quantities are represented 



by partition functions in the forward case and distribution amplitudes in the case of non-forward 
scattering. Here, and for the description of the evolution equations below, the general form of 
the operators O t (kix, k^x) has to be taken into account. It will be shown that their Fourier 
transforms with respect to K_xp + and define the above mentioned distribution ampli- 

tudes. For the following considerations it is also useful, as it was in the twist decomposition, to 
study the non-local operators temporarily to arbitrary values of x and to perform the light cone 
limit only afterwards. 

In the following we consider the contributions of twist-2 only. In the case of purely electro- 
magnetic interactions the matrix elements at lowest order in the coupling constant, using the 
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equations of motion, can be represented by two basic hadronic matrix elements, including also 
the spin dependence: 

Q D (x;p 2 , Pl ) = x»Q°(p 2 ,p 1 )=x»u(p 2 ,S 2 )j IJ ,u(p 1 ,S 1 ) (4-1) 
Q P (x;p 2 , Pl ) = x^Q p (p 2 ,p 1 ) = x^u(p 2l S 2 )j I a llu p^u(PuSi) , (4.2) 



with 



i 



M is the nucleon mass, and u(p, S) denotes a free hadronic Dirac spinor. Here the indices D 
and P refer to the Dirac and Pauli structure, respectively. 

Let us first consider the matrix elements of the scalar twist-2 quark operator. Because of 
translation invariance, 

(p 2 \O r (K 1 x,n 2 x)\p 1 ) = (p 2 \O r ((n + - k_)x, (k+ + K-)x)\pi) 

= e iK+xp -(p 2 \O r (-K„x,K„x)\ Pl ) , (4.3) 

it is completely sufficient to consider the matrix elements of the centered operators: 

(p 2 |O a '* wiBt2 (-K_a:,K_a:)|pi) = ifj(K_xp + , K 2 _x 2 , Pl p 2 , fi 2 )Q J (x;p 2 , Pl ), (4.4) 

where the summation is taken over J = (D, P). If considered independently these functions will 
be denoted by fo = g and fp = h. The generalization to arbitrary values of Ki is trivial due 
to translation invariance, Eq.( |4.3j ), through which a factor exp{in + xp-} occurs on both sides of 
Eq. ( |4.4D . The projection onto the light-cone immediately gives 

( P2 |09' twiBt2 (-K_£,K_z)|p 1 ) =if q J (K_5;p + ,K„xp„,p 1 p 2 ^ 2 )Q J {x]p 2 ,p 1 ) . (4.5) 



The twist-2 vector operators, because of the additional derivative in Eq. ( |3.8| ), do not have 
such a simple representation. Using that relation between the vector operator and the scalar 
operator for general arguments we obtain 

(p 2 |0«' twtot2 (K 1 £,K 2 5)|p 1 )= f 1 d\d^^ 2 \0^ t2 {^x^ 2 \x)\ Pl ) (4.6) 



= i J d\d* (fj(K^\xp + ,K„\xp^,K 2 _\ 2 x 2 ,pip 2 ,/j, 2 )Q J (x;p 2 ,pi)e iK+Xxp - 

Let us now perform the Fourier transformations of fj with respect to K-\xp± = K,xp±, 

f(Kxp + ,Kxp_,K 2 x 2 ) = J Dze- iKx{p * z * ) f{z + ,z„,K 2 x 2 ) , (4.7) 

where 

{piZi) = p(z) = p + z + + p_z_ with z± = \ (z 2 ± zi) , (4.8) 

defines a scalar product in the space of two-vectors labeled either by (1, 2) or (+, — ). Also these 
functions, if projected onto the light-cone, are entire analytic functions with respect to n^xpi 
and therefore the support of their Fourier transforms is restricted to the range 

- 1< z, < 1 . 
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Thereby, we assumed that the kinematic decomposition of the matrix elements does not introduce 
additional kinematic singularities. Therefore, the integration measure in Eq. ( |4.7| ), analogous to 
Eq. ( |2.13| ), is given by 



Dz = -dz 1 dz 2 9(l-z 1 )9(l + z 1 )9(l-z 2 )9(l + z 2 ) 

= dz + dz-9(l + z + + z-)9(l + z + - zJ)9{\ - z + + zJ)9{\ - z + - z_) 

The A-integration in Eq. (^4.6| ) can be performed leading to 

(p 2 \Oy wist2 { Kl x, K 2 x)\ Pl ) = i J Dze- iK -^ )z * 



(4.9) 



(4.10) 



x {Fj(z + , z_, K + xp„, k 2 _x 2 ) (Q^{p2,Pi) ~ iK- P »(z)Q J '(x;p 2 ,p! 

+ (iK + p^.d iK+xp _ +2x ll K 2 _d K 2_ x 2^ F q j(z + , z-, K + xp-, k 2 _x 2 )Q j (x; p 2 , pi)} 

Here, the following functions F q (z + , K + xp^, k?_ x 2 ) are introduced: 

, 2 2 \ f 1 d\ n f z + Z- 

r [Z + , Z-, K + xp_, K_X ) 

F (-£-!-, ? fc-\-XP— , K X ) 



A 2 



A ' A 



-/ -± ^, k_X x 6(A - \z + \)G(\ - \z-\)e 



0, 



IK + Xp- 



F) (z + , Z-, K + Xp^, K 2 _X 2 ^ 



(4.11) 
(4.12) 



1 ^/ (£+ £z K 2_ AV j Q(A _ | 2+ |)e(A - |*_|)e 



A 



A ' A 



(S^ ^f) (* + , k 2 __x 2 ) (4.13) 

£dXd KW f (f±,^, K 2 _A 2 x 2 ) 9(A- |2+|)6(A- |z_|)e^- . 



As will be shown in Section [5] one consequence of the representation ( |4. 1 1|) are the integral 
relations in polarized deeply inelastic scattering in the forward case |B| and [2D], as well as the 
Callan-Gross relation [pTJ . For p_ = one obtains: 



1 dz 



Thereby the derivation f'(z+, Z-, k 2 _x 2 ) = (d K 2 x 2f) (z + , 2_, k 2 _x 2 ) was introduced in Eq. (|4.13| ). 
Furthermore, a third Fourier transformation can be performed either with respect to n 2 _x 2 |TT| 
or, as will be used here, with respect to which allows to study the (l/Q 2 )-dependence of 
the distribution amplitudes 



f(^xp + ,K^xp„,K z _x z ) = J Dz J d*qe- iK - {xp ^e- iK - qx f(z + ,Z-,q 2 ) . 

In the limiting case of forward scattering we have: 
Q D (x;p,p) = u(p)x M 7 M w(p) = 2xp and Qp(x;p,p) = u{p 2 )jjx^a IJ , l/ p l Lu{pi] 

and one obtains (g = fr>) 

(j,\O q > twisa (-K-x,K-x)\p) = i2(xp) g q {2K_xp) 



(4.14) 



(4.15) 



i2(xp) J dz+e- 2iK -^P>+ 



+1-1*4-1 
-l+|2+l 



dz_g q (z + ,z_) . 
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Analogous decompositions and representations are valid for the matrix elements of the pseudo- 
scalar and pseudo-vector twist-2 quark operators. The corresponding partition functions are 
denoted by fj 5 = (gf, hi) and Fj 5 = iff), respectively, and the Dirac and Pauli structures 
are to be replaced by , Q§ instead of Eqs. ( |4. 1| , 4.2 ) , substituting 'y jU — ► 'YsTV 1 

and a^ v — > 750^, 

respectively. 

Let us now extend our consideration to the twist-2 scalar and vector gluon operators, 
Eqs. ( |3.14j ) and (|3.15| ). Completely analogous to Eq. ( |4.4|) the matrix elements of Eq. ( |3.14|) 
are written as 

{p 2 \0 G ' tmst2 (-K_x, K_x)\pi) = i(ixp + ) fj W (K^xp + ,K_xp^ 1 K 2 ix 2 ,piP2,iJ 2 )Q J (x;p2,Pi) ■ (4.16) 

The additional factor ixp + has been introduced in order to compensate for the different K-scale 
dimensions, d q = 1 and do = 2, of the quark and gluon operators, respectively, cf. also Section 
6.1. 

The matrix elements of the twist-2 gluonic vector operators, Eq. ( |3.15|) , are obtained in the 
same manner as for the corresponding quark operators: 

(p 2 \ G >^( Kl X,K 2 x)\ Pl ) = i J Dze -i K -(* Pi )zi (4.17) 

x \{ixp + ) Ff {1 \z + , z_, K + xp-, k 2 _x 2 ) (Q^(p2,Pi) - iK^p fl (z)Q J (x;p 2 ,Pi) 
+ (iK + P-d iK+xp ^ + 2x fM K 2 _d K 2_ x 2^j Ff (1) (z + , z-, K + xp_, K 2 _x 2 )Q J (x;p 2 ,p 1 ) 
+ ip^Ff {1 \z + ,z_, n + xp-, K 2 _x 2 )Q J (x;p 2 ,pi)\ _. 

J x=x 

We remark that in our preceding paper, Ref. ||13|| , for simplicity we introduced formal defi- 
nitions for the matrix elements of the scalar quark and gluon operators which are sufficient for 
the derivation of evolution equations. There we introduced the following two representations 
referring to different kinematic variables, 



(p 2 |0^ ist2 (^ 2 )bi> = e iK +* p -(ixp + ) j Dze- iK -^ Zi f^(z + ,z.), (4.18) 
{p 2 \0 G ^ t2 (K hK2 )\ Pl ) = e iK ^-{ixp + f J Dze- iK -^ Zi f^(z + ,z^), (4.19) 

which are now to be replaced by the explicitly determined expressions Eqs. (|4.4j ) and ( [4.16|) 
projected onto the light-cone. In the limit of forward scattering these functions describe the 
hadronic parton densities. 

However, contrary to the case of the quark operators, for the matrix element of the gluon 
operator there exists no natural kinematic decomposition. A second parametrization used earlier 
is 



k_ (p 2 \O q ' twist2 (K 1 ,K 2 )\p 1 ) = e iK+ip - J Dze- iK - ( * p ' )z >tf 2) (z + ,z_), (4.20) 
k 2 _ {p 2 \0 G ' twiBt2 {K 1 ,K 2 )\pi) = e iK+ip - J Dze- iK - (ip ' )z >ff 2) (z + ,z_) (4.21) 
corresponding to 

«- (P2\0 q ' twmt2 (-K^x, K-X)\px) = (xp + )~ 1 ff 2 \K^Xp + ,K_Xp^,p l p 2 ,jJ, 2 )Q J (x;p 2 ,Pi), 
(p 2 \0 G ' tW1St2 (-K-X, K-X)\pi) = (xp + )~ 1 f ( f {2 \K-Xp + , P1P2, fJ 2 ) Q J {x] p 2 , Pi) ■ 

(4.22) 



13 



Because the vector operators are not used in the following we do not write down the respective 
representations for their matrix elements. 

Whereas the first choice in the case of forward scattering leads to the correct matrix elements 
being proportional to (2xp) 2 , the second choice leads to a much simpler form of the anomalous 
dimensions. Note, that in both cases, up to a trivial reduction of powers of k_ on both sides, 
everywhere the combination k_x appears as a natural variable. 

Let us point to the fact that for both the twist-2 quark and gluon operators the (pseudo)vector 
distribution amplitude Fj, Fj$ are directly related through Eqs. ( |4.11| ) and ( j4.12|) to the 
(pseudo) scalar distribution amplitude fj, fj$. Therefore we need not to investigate the renor- 
malization properties of the (axial)vector operators independently: It is sufficient to work out 
the renormalization of the (pseudo) scalar operators and to derive the renormalization group 
equations for them only. 



5 Compton scattering amplitude in leading order 

With the prerequisites provided in the preceding Sections we now derive the asymptotic repre- 
sentation of the Compton scattering amplitude. Again, we consider the contributions of twist-2 
operators only and we restrict the analysis to leading order. Starting from the matrix elements 
of the non-local operators which depend on the two independent variables z± we arrive, af- 
ter performing the Fourier transformation of Eq. ( |1.2| ), at double-variable distributions. 4 In 
the literature also single-variable non-forward evolution equations are studied by imposing an 
additional kinematic constraint. The latter description is suited for special kinematic situations. 



5.1 Double- variable distributions 
5.1.1 Non— forward scattering 

To obtain the double-variable distributions we use the decomposition Eq. (|4.10|) of the Fourier 
transform of the Compton amplitude 

T llu (p + ,p_,q)=i J d*xe iqx (p 2 ,S 2 \RT(J ll (x/2)J l/ (-x/2)S)\p u S 1 }. (5.1) 



Afterwards we have to carry out the limit x — > x. In the r.h.s. of Eq. ( |5.1|) the representation 
of Eqs. ( |2.12| - |2"TT7| ) is used. Furthermore the matrix elements are expressed in terms of hadron 
states, Eqs. ( [Q| , f4.6|) . The arguments k± of the coefficient functions C a ( 5 )(x 2 , k + , k_, /i 2 ) in 
Eqs. (2.14, |2~T5| ) are fixed at k + = and = ±1/2. For the representation of the Compton 



amplitude we could use the matrix elements at these points. However, for the investigation of the 
evolution of the matrix elements their representation at general values of k + and K- is needed. 



The Fourier transforms of the individual terms are obtained by, cf. [^5] 



e g 



/ 



27r 2 (x 2 — is) 2 q 2 + is 

d 4 x ■ x v x" g vix q u q^ 

e q ^-n i To = -i n . +2i 



2tt 2 (x 2 — is) 2 q 2 + is (q 2 + is) 2 

4 The general Lorentz-structure of the Compton amplitude was investigated in Refs. ^H, cf. also [p3| . 
Real-photon processes were considered in Ref. |24] before. In general the Compton amplitude consists of 18 basis 
elements. For qf — q\ four spin-independent and eight spin-dependent distribution amplitudes contribute ||. 
Their number reduces to two and four, respectively, for q\ = q\ =0. 
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d 4 x 



,iqx_ 



2tt 2 (x 2 - ie) 2 (q 2 + ie) 2 (q 2 + ie) 3 ' 

In the generalized Bjorken region one obtains, noting that £ = Q 2 /qp + , rj = qp-/qp + , 



T;r t2 ( P+ ,p-,q) 



Dz 
+ 



£ + 2 + + 7],2„ — £ — z + — t]z_ — ie 
1 1 



(5.2) 



+ 
+ 



(f + 2+ + ^- 
1 


— ie' 




(£- 


^+ - 


T}Z- — 
1 


ie) 2 ) 


(f + z + + rjz- 
1 


— ie] 


3 


(£- 


^+ - 
1 


T]Z- — 


ie) 3 ) 


£ + Z + + T]Z_ 
1 


- ie 


V 




— 772:. 


. — ie) 
1 




— ie] 


2 




z+ - 


— 


ie) 2 ) 



,(2) 



In this representation of 7)^ we refer to the essential terms only. One may generalize Eq. 
accounting for current conservation explicitly. For the case of forward scattering this represen- 
tation was given in Ref. p6|. 



The functions (p + ,p-, q; z + , z_) being (anti) symmetric w.r.t. the Lorentz indices are 



given by 



fii/ 



771(1) 
5,fiu 



r (2) 



qp^ 



Qt(P2,Pi)Fj(z + ,z_), 



(5.3) 



/ 1 



2 r 



9^qp{z) - {q^p v {z) + q u Ptt{z)) 



q a Qt(p2iPi)F>j(z + ,z_), 



qp A 



\qp. 

ie»u aP q a —Q J /F J>5 {z + ,z. 

qp+ p 



q a Qt(P2,Pi)Fj(z + ,z_ 



-ie^ q^ppiz) 



qp^ 



q a Q J a 5 ( P 2,Pi)Fj, 5 (z + ,Z- 



(5.4) 
(5.5) 
(5.6) 



where p(z) is defined in Eq. (|4.8|) . Let us remind that the functions i*j,(5) are taken in the limit 
x 2 — > and that Fj(z + , z_) results from the Fourier transformation of d K 2 X 2F. The distribution 
amplitudes -Fj,(5) and Fj are complex quantities in general. The representation of the LCE, 
Eq. (|2.3|) , implies due to its construction that the current conservation for the incoming and 
outgoing photon 



q 2 ,T^ = 0, T^q lu = 
is not obeyed in explicit form yet. However, one may impose these conditions. 



(5.7) 



5.1.2 The forward scattering limit 

In the limit of forward scattering, p_ = p 2 — p± — > 0, rj = 0, the Compton amplitude does not 
depend on the distribution amplitudes Fp t (5)(z+,zJ) anymore. Using the normalizations 

u(p)u(p) = 2M, u(phtiu{p) = 2p^ 
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for the free hadronic Dirac spinors we obtain for the symmetric part of the Compton amplitude 

1 



T-isym 



9 



pq 



£ + z+-ie z + -ie 
1 



+ 



1%(lv J Dz 



(£ + z+-i£) 3 (£-z+-i£) 3 / 
The ^.-integral can be performed 



G' q (z + , z 



(5. 



f+l-\z+\ /-+J--|*+I /• 

/ = / cfe_ / 

J-l+lz+l J-l+\z+\ JO 



2-|- li 



(5.9) 



After partial integration of Eq. (|5.8|) , 

/•Sign(z + ) ^ r+l l [ /-sign(z + ) ^ 



-i 



(£±^ + -^) 2 
we finally obtain 



T-isym 



1 



1 



^ + z+ — ie £ — z + — is _ 



- ze) 3 (f - z+ - is) 



(5.10) 



In the same way we can treat the antisymmetric part of the Compton amplitude. Using the 
normalization condition 



uiphslpuip) = -2S>, S 2 



-M 



we obtain analogously 



T-iantisyn 



IS 



fin 



7/3 W 

(pqf 



+1 



—IE 



IP 



q,S y+i 



(pg) 7-1 



1 1 

+ 



£ + z + — is ^ — z + — is 



sign(z+) dz 

2-|- Z 

sign(2+) ^ 



/ — 0f(*) 



(5.11) 



5.1.3 Relations between structure functions 

The hadronic tensor is related to the Compton amplitude for forward scattering by 

Wfa, = — ImT^ . 



(5.12) 
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In the case of purely electromagnetic interactions the Lorentz-structure of the hadronic tensor 
is 



w, 



„, = -g^ix, Q 2 ) + ^f 2 (x, g 2 ) + ^f 4 (x, g 2 ) + + 

pq pq pq 

Q X S a , rfl ... q X (pqS* - SqjT) 2 



F 5 (x,Q 2 



(5.13) 



where x = xb] = £■ Whereas the polarized structure functions in Eq. ( |5.13| ) form already the 
minimal set, current conservation 



q»W, u = W^q v = 

relates the structure functions F4 and F 5 to F\ and F 2 by 

1 



F 4 (x,Q 2 ) 
F 5 (x,Q 2 ) 



Ax 2 
1 



2xF 1 (x,Q 2 ) + F 2 (x,Q 2 ) 



2x 



F 2 (x,Q 2 ) . 



(5.14) 

(5.15) 
(5.16) 



The forward Compton amplitude Eqs. (|5.8|) and ([5 -11 ) are related to the hadronic tensor 
Eq. dU3|) using || 



lim 



p-L±(-i)»_^ s^h 

x n v ' (n - 1)! 



x 



(5.17) 



e ^o+ (x ± iey 

We now assume that g? 5 Jz+) are real functions. The tensor structure of Eq. ( p.lQ) implies 

F 1 (x,Q 2 ) =F 5 (x,Q 2 ) (5.18) 
at leading order, from which, cf. Eq. ( {5.16j ), the Callan-Gross relation pl | 

F 2 {x,Q 2 ) = 2xF l {x,Q 2 ) (5.19) 

follows. Furthermore, one finds 



F,(x,Q 2 ) = -F 1 (x,Q 
x 



(5.20) 



which relates the distributions G' q (z+, z_) and G q (z + , zJ) for forward scattering at leading order. 
The structure function F\ obtains the representation 



F 1 (x,Cf) = -\§<(x)-F(-x)] 



(5.21) 



We suppressed, as in the preceding Section, the scale dependence of the functions g%-, . In lowest 
order of perturbation theory the g 2 -dependence of the structure functions results from the 
identification fi 2 = Q 2 . One may identify this representation with results being obtained in the 
quark-parton model at leading order by defining 



Nr 



g°(x) 



e ^O) 

9=1 

N f 



q=l 



(5.22) 
(5.23) 
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Here q(x) and q(x) denote the unpolarized quark and antiquark densities, respectively, and e q is 
the electric charge of the quark. 

We now turn to the polarized case. The Lorentz structure of the two tensors in Eq. (|5.11 ) 
refers to the structure functions g\+g 2 and g 2 , respectively. After performing the ^-integration 
one obtains for the twist-2 terms the Wandzura-Wilczek relation |19| 



g 2 (x,Q 2 ) = -g^x.Q 2 ) + / — g x (z, Q : 

J x Z 

gi(^Q 2 ) = -J2e 2 g [Aq(x) + Aq(x)} 

A 9=1 



where 



and we identified 



05 0*0 



E 

9=1 



e q Aq(x) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



9=1 



In the case of general electro-weak couplings also a second integral relation, Ref . |21J , is obtained. 
For both these relations the parameter-integrals relating the vector and scalar operators which 
are discussed in Section ^ form the background in the present approach. The new derivation of 
these relations outlined above could circumvent their usual derivation by means of the moments 
in the local operator product expansion and leads to the integral form directly. We note that 
recently also three integral relations for the twist-3 contributions to the polarized structure 
functions 9j\j = i for the case of forward scattering were found in lowest order in the coupling 
constant, cf. Ref. [26|. Unlike the case for twist-2 
requiring an all-order resummation of the (M 2 /Q 2 )-terms. 



20] all the relations are integral relations 



5.2 Single— variable distributions 

Up to now we have considered representations which contain two distribution parameters z + , z-. 
In some connections another representation may be helpful. As an example let us introduce in 
representation Eq. ( |5.2|) a new variable t by z + = t — T]Z-. Then the integrations factorize and 
we can write 

zrw-,g) - fyiisrh-e - izh*)^*-'**^ (5 - 28) 

\{t + t-ie) 2 ^(S-t-ie) 2 ) » v \(C + t-ief (Z-t-ie) 3 ) ^ 

, ( 1 , pd) , ( 1 

\Z + t-ie t-t-ie) \(£ + t-ie) 2 (£-t-ie) 2 ) 5,/Xi/ 

Also within the functions F^) (p + ,p_, q, Fj(z + , z^;pip 2 , /i 2 ); /x 2 ) the variables z + , z_ have to be 
transformed thereby changing these functions to Fjrg\(z+ = £_ — r)Z-,Z-.). Then the remaining 
^.-integration can be performed to give 

r+l 

Fj(5)(t,v) = J dz-G{l + t-riz- + z-)e(l+t-riz-- zJ) 

x6(l - t + rjz~ + *_)e(l - t + rjz- - z^)F J(5) (z + = t - r]z_, zJ). (5.29) 



18 



These new partition functions Fj and Fj^ depend on the variables t and 77 and therefore we have 
F^lfiviP+iP-iTitiV^ 2 ) = F (5)^u (p+^P-^^J^V^ 2 )^ 2 )- Obviously, in the limit of forward 
scattering we obtain F^ ^(2p, 0, q; t, 0; /i 2 ). Note that the partition functions F, F' and F 5 are 
directly related to the partition functions /, /' and f$ which are defined by the matrix elements 
of scalar and pseudo-scalar twist-2 quark operators, Eq. (|4.4|). This demonstrates again that 
for the description of the virtual Compton scattering, at least in leading order, the properties of 
the scalar operators are sufficient. 

The renormalization of the twist-2 operators will be considered in the following Section. The 
evolution equations of their expectation values are dealt with in Section [7]. 

6 Operators of twist— 2 and their anomalous dimensions 

In this section we discuss the renormalization properties of the twist-2 light-ray operators. As 
has been shown above we can restrict these studies to the case of scalar operators for which any 
free tensorial index of the various operators is saturated by multiplying with as many as necessary 
vectors of x^. Moreover, there holds the much wider statement that the operators 0" ,twlst2 and 
O a,twlst2 are completely equivalent concerning their renormalization properties. This can be seen 
in the following way. We consider the renormalization group equation of the vector operator, 

/^ 2 A^' twist2 (KiX,K 2 x;/i 2 ) = r i ^W47(/t 1 ,/t 2 ;4,4;/i 2 )0^ twist2 (4x,4x;/i 2 ) , (6.1) 
d[i l Jk 2 

cf. Eq. ( |6.20D . Obviously, the scalar operator 

O a ' twist2 («i, k 2 ) = £ ,7 0^ twist2 (/« 1 , k 2 ) (6.2) 

fulfills exactly the same renormalization group equation with the same anomalous dimension 
kernel 7(/Ci, k 2 ; k' 2 ), since the multiplication of both sides by x u commutes both with the 
differentiation on the left and with the integration on the right hand side. The same conclusion 
can be drawn for the gluon operators. 

Because the renormalization properties of the various twist-2 operators are independent of 
their special Lorentz structure we restrict ourselves to the scalar operators. Contrary to the 
previous considerations we will account for the flavor content of the operators because in the 
flavor singlet case twist-2 quark and gluon operators mix under renormalization. To simplify 
the subsequent considerations, and the necessary computations of Feynman diagrams, we apply 
the axial gauge where the phase factors are U(k\, k 2 ) = 1. 

We consider the following scalar flavor non-singlet (NS) and singlet light-ray operators 5 : 

O ns (ki,k 2 ) = x^ ^ ^(n^Xf^^^x) - ^(K 2 x)\ f ^^(Kix) (6.3) 

0^ s (k 1 ,k 2 ) = x^ \ ^(^1^)^/757^(^2^) + i){K 2 x)\ f ^^^(Kix) (6.4) 

O q (K 1 ,K 2 ) = x^ ^ ^(kix)^^1)(k 2 x) - ^(^X^^KxX) (6.5) 
G ( Kl , k 2 ) = x»x v \ [F^K^F^x) + F^i^F^x)] (6.6) 



J Possible trace terms vanish and the general dependence RT[0]S on the renormalization procedure is under- 
stood but has been omitted here. 
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Of (ki, k 2 ) 



^ 1 



(6.7) 
(6.8) 



Here A/ denotes the generators of the flavor group SU(Nf), where Nj is the number of active 
quark flavors. The scalar operators Eqs. (|6.3| , |6.5|) and (|6.6| ) contribute in the case of unpolarized 
Compton scattering, whereas the pseudo-scalar operators Eqs. ( |6.4| , |677| ) and ( |6.8| ) are relevant 
for Compton scattering off polarized hadrons. 



6.1 General properties of non-local anomalous dimensions 

We consider now some general properties of the anomalous dimensions of these operators which 
hold at any order of perturbation theory and are relevant in the following. In order to cover the 
general case we refer to the scalar singlet operators denoted by A with A = (q, G). R 
The renormalization group equation under consideration reads: 

li 2 -—O a (kix,k 2 x-,li 2 ) = ( dK[dK 2 ^ AB (Ki, k 2 , k[, k 2 ] li 2 )O b (k[x, k' 2 x; li 2 ) . (6.9) 
a fx J K2 

7 AB denotes the non-local anomalous dimension matrix which, through the strong coupling 
constant a s (/i 2 ) = g 2 s (li 2 ) / (4ti) , depends on the renormalization scale Analogous relations 
hold for the pseudo-scalar case. The non-singlet cases are covered by the respective projection 
in flavor -space of the quark-quark submatrix. 

Since the variables K{ simply parametrize points on the light-ray these anomalous dimensions 
are invariant under reparametrizations p8|, i.e. under translations and scale transformations,^ 



7 AB (ki, k 2 ; k'u k 2 ) = 7 AB (/t 1 — k, k 2 — k; k' x — k, k 2 — k) (6.10) 

= A dAS 7 AB (Afi: 1 ,A/t 2 ;A4,Afi: , 2 ) , (6.11) 

where 

d AB = 2 + d A -d B} (6.12) 

with d A being the difference of the canonical and the /t-scale dimension of the operators A . 
The latter is given by the number of factors x emerging in the operator A , i.e., 

d q — 1 and da = 2 . (6.13) 

Therefore, the number of relevant ^-parameters may be reduced by two. We consider two 
(related) choices in order to normalize the anomalous dimension matrix, shifting by either K\ or 
n + and scaling by either (k 2 — Ki) _1 or (k_) _1 , respectively. One obtains 

(k 2 - « 2 ; /<, 4) = 7 AB (0,l;ai,l-a 2 ) = i^ B («i,« 2 ), (6.14) 

4(«_) <,Afl 7 AB («i,« 2 ;«i,/4) = ^ AB (-l,+l; Wl ,w 2 ) = K AB (w 1 -w 2 ,w 1 + w 2 lQ.15) 



6 The respective relations hold synonymously for the three possible cases of non-singlet evolution equations j27j 
which have to be suitably projected out. They are all equivalent in leading order, and two of them are even 
equivalent in next-to-leading order. 

7 In the remaining part of this and the following section the explicit dependence on x and /i 2 will be suppressed 
in the operators, anomalous dimensions and evolution kernels. 

8 These properties of the anomalous dimensions may be traced back to the independence of the renormalization 
properties of the operators from shifting or scaling their arguments. 
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where K i Kl K 2 K 2 (a 1C \ 

w c c «i = , —a 2 = , (6.16) 

k 2 — Ki k 2 — K\ 



k' + — /«_■_ K 1 



W\ = «l — «2 = , W2 = 1 — «1 — «2 = , (6-17) 



respectively. The variables k' are related to the parameters a and w by 

/VA //ti(l - «i) + /t 2 ai\ 1 //«l(l — lOl + W 2 ) + K 2 (l + Wi — W 2 ) 

\K 2 J \K 2 (l - Ol 2 ) + Kia 2 J 2\K 1 (l—W 1 -W 2 ) + K 2 (l+Wi + W 2 ) 



(6.18) 




Figure 1: Integration ranges of the parameters (ai,a 2 ) and (wi,w 2 ), Eqs. ( |6.21| , |6.22|) . 



Therefore, it seems to be quite natural to introduce instead of the variables (k[, k' 2 ) as new 
variables either (ai,a 2 ) or (wi,w 2 ). The integration measures are related by 



Dk' 



[k 2 - KxfDa = A(k_) 2 Dw , 



where 



Dk' 

Da = D a + D x a 



«2J 



dKidK 2 0(/«i — K^)G(/Cj — K2)0(^l — K 2)©( K 2 

daida2 0(ai)0(a 2 )0(l — ot\ — a 2 ) 

+ 6(1 -ai)e(l -a 2 )0(ai + a 2 - 1)1, 
Dw = D + w + D^w = | dwidw 2 0(1 + u>i — w 2 )Q(l — W\ — w 2 )Q(w 2 ) 

+ 0(1 +wi +w 2 )Q(l - wi + w 2 )Q(-w 2 ) 



(6.19) 

(6.20) 
(6.21) 

(6.22) 



The w-representation was derived in Ref. || for all orders in the coupling constant. The measures 
Da and Dw have been divided into two parts related by <-> 1 — and w 2 <-> -w 2 , respectively. 
Usually only the first part (unshadowed in Fig. 1) is considered in the literature. Because of 
symmetry requirements for the anomalous dimensions being discussed below this more general 
representation proves, however, to be appropriate. 
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Let us now consider the symmetry of the quark and gluon operators, and likewise of their 
anomalous dimensions, under changing (k x , k 2 ) ( K 2; K i) and k' 2 ) <-» (k 2 ,Ki). Looking at 



Eqs. (|6.5| -|6.8| ) it is easily seen that 

O a (k u k 2 ) = (-1)^O a (k 2 , Ki ), (6.23) 
Oi^K,) = -{-l) dA A {K 2 ,K X ), (6.24) 

holds. We first consider the scalar case. From Eq. (|6.23|) we obtain 

j AB (K hK2 ; k' x ,k' 2 ) = (-l) d «7 Afl («2,/ci;/4,«i) = (-l) d Wi,^; 4, «i) • (6-25) 
Using reparametrization invariance, Eqs. ( |6.10|) and ( 6.11 ), we obtain 

K AB (a u a 2 ) = ^ AB (a 2 , ai ) = (-l) dB K AB (l - a u 1 - « 2 ), (6.26) 
K AB (u?i,u? 2 ) = K AB (- Wl ,«; 2 ) = (-l) fe ^ B ( Wl ,-«; 2 ). (6.27) 

The renormalization group equation (|6.9| ) transforms into 

H 2 -^O a (k u k 2 ) = j Da{n 2 -n l ) dB - dA K AB {a ll a 2 )0 B {K' l ,^) 1 (6.28) 



^^0 A ( Kl ,K 2 ) = J Dw{K_) d -- d -K AB {w l ,W 2 )0 B {K' l ,K l 2 ) 



1 1— W2 



! 



dw 2 J d«; 1 («_) fe - d ^^(«; 1 ,«; 2 )0 B («i,4), (6.29) 



-1+1U2 



using the kernels ^^(ai, a 2 ) and ^^(iwi, w 2 ), Eqs. ( |6.14|, |6.15|) .P1 In the second line of 
Eq. ( |6.29|) we have taken into account the above decomposition ( |6.22| ) of the integration mea- 
sure and the second of the equalities ( |6.27| ) to restrict the integration onto the usual range by 
symmetrizing the anomalous dimension kernels in an appropriate way: 

K A l{w u w 2 ) = i[K AB ( Wl , W2 ) + (-l)^^( Wl ,- W2 )] . (6.30) 

Here K AB (wx, w 2 ) is defined in the range of the measure D + w only. Even more, with the help of 
the first of the equalities ( |6.27|) the -u^-integration can be restricted to the range < W\ < l — w 2 . 
Using these kernels we can write the renormalization group equation for the centered scalar 
operators in terms of the variables w as follows: 

1 1—W2 

Ii 2 — 2 0\-k„,k„) = 2 Jdw 2 y'diUiKi fl -^^(tw 1 ,iu 2 )0 B ((«;i-«; 2 )«_,(«; 1 + «; 2 )K_). (6.31) 

" 

In terms of the a-variables these expressions are more lengthy. 

In the pseudo-scalar case the sign in the last expression of Eqs. ( |6.25 - 6.27] ) , as well as in 
Eq. (|6.30|) has to be changed into — (— l) dB in accordance with Eq. (|6.24j) , which also implies the 
behavior of the sign functions a^, defined in Eq. ( |6.52|) below. The corresponding kernels will 
be denoted by A7, AK and AK, respectively. 

9 To get rid of the unwanted K-dependent factors in these equations, it seems to be reasonable to multiply 
the scalar quark and gluon operators O a by {k 2 — Ki) dA or k_ if the a- or ^-representation has been chosen. 
However, this does not really resolve the problem since an additional factor (1 — a\ — at2)~ dB or u>^ ds is introduced 
which multiplies the kernels K and K, respectively. 
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6.2 Relations between non— local and local anomalous dimensions 



Obviously, the general statements proven here for the non-local operators and their anomalous 
dimension kernels can be formulated for the local operators and their anomalous dimension 
matrices, too. With the definition Eq. (|2.9|) of the local operators 0^ in2 the renormalization 
group equations read 



1 " qA 



2 ~ n\ri2 



E AB qB 
mi,n.2i n',,ni 2 n' 1 n' 2 



(6.32) 



i 1 , n 2 



with 



4B 



Qrn Qn 2 i- Kl 



I ( K l) Ul {^Y 2 AB ( I i \ 

7 (ki,k 2 ;ki,k 2 ) 



n'A 



K1=K2=0 



IK 2 I Op 112- 

First of all we remark that the symmetry relation ( |6.25|) is simply transformed into 



AS 
/ 711,712; n^ni^ 



' _\ \d AB AB 

\ J In2,ni;n' 2 ,n' 1 



d B AB 

'711,712; n 9 ,n. 



(6.33) 



(6.34) 



Furthermore, taking into account the reparametrization invariance, Eqs. ( |6.14|) or (|6.15|) , together 
with the corresponding representation of k[, it is easily seen that the values of and n\ are 
restricted by the relations 



N = rii + n 2 = n\ + n' 2 + d B - oIa 



(6.35) 



since otherwise the anomalous dimension vanishes. From this it is obvious that only those local 
operators mix under renormalization which are specified by N. Therefore, we may choose 



rii = N — n, n 2 = n; n[ = N — n' , n' 2 = n' — d B + d>A , 
and the local anomalous dimension can be written as 



(6.36) 



AB 



Tni, 



712 ; ,n' 2 



AB 



lN-n,n;N-n' ,n' -d B +dA 



N AB 
i n n' 



(6.37) 



Now we go over to the variables k± in place of (k\, /t 2 ) which are favored because of the relations 
( 6.15| ) and (|6.17 ). For simplicity we perform the next steps for the non-singlet local anomalous 
dimensions only. Analogous results are obtained in the singlet case. Let us define 



d n J 0(«ix,« 2 x)| K 



711,712; n^n^ 



no 



±=0 ' 
/ \n'. I .J \ni 



Dw 



Ho 



K(w 1 ,w 2 ] 



K±=0 



(6.38) 
(6.39) 



with the abbreviation d± = d/dn±. Taking into account the representation of k' ± by Eq. ( |6.17|) 
and decomposing (k' + ) Wi = (k+ + WiK.)" 1 we obtain 



N-n an N—n' 



N 



^nn 1 — TVi 



QN-n Q 



(J- \- (N-n'\ 
a,„/| 2^ { l ) K 



N-n'\ N-n'-lJ+n' 



(N — n')\ n'\ h n 
(£) f Dwwr n 'w^K( Wl ,w 2 ) . 



Dww{w% K(w 1 ,w 2 ) 



(6.40) 



From this equation some remarkable facts can be read off which hold at any order of perturbation 
theory. First of all, the local matrix of anomalous dimensions is triangular, i.e. n' < n and, 
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secondly, it is universal, i.e., independent of N. Thirdly, because in the case of forward scattering, 
p_ = 0, the variable k + disappears from the expressions, cf. Eq. ( |4.3| ), the relation N = n = n' 
is implied with the consequence that the diagonal elements of that matrix are the forward 
anomalous dimensions. In the unpolarized case one obtains 

7n = |(l - (-1)") f Q dw 2 w r i 1 2 dwiK (wi, w 2 ) j . (6.41) 

The anomalous dimensions vanish for even values of n.Q 

Since the local anomalous dimensions in the forward case are moments of the splitting func- 
tions, the latter are given by 

P(z) = 2 f 1 * dw x K Q {w u z) (6.42) 



performing the inverse Mellin transformation. Here z denotes the partonic momentum fraction. 
Analogous relations hold in the pseudo-scalar case, with the difference that now the even values of 
n are to be taken. These considerations show the value of the w-representation, which naturally 
introduces a partial diagonalization in the variable w 2 through the Mellin transformation. It 
becomes also evident that the Wi-dependence is related to the off-diagonality of the anomalous 
dimension matrix. 

Let us now consider the local anomalous dimensions for the singlet case. The local scalar 
light-ray operators are given by 

N 6 q n = d^d^Urx,^) with 0<n<N, (6.43) 

K±=0 

N G _ l = ^- n 9r 1 G (/t 1 x,/t 2 x)| K± _ o with l<n<N, (6.44) 

where N + 1 equals the total spin of the local operators mixing under renormalization. The 
difference between quark and gluon operators results from the fact that their K-scale dimensions 
equals d q = l,da = 2. Obviously, n Oq has no counterpart to mix with. This situation can be 
resolved if 

(Q q ,0 G ) = (<9_0 9 ,0 G ) (6.45) 
is used instead of (O q ,0 G ). The corresponding renormalization group equations read 

^/Ot, = E N lt^ n ,-i N Q%-x with l<n,n'<JV, (6.46) 

a f i n'=l 

where the local anomalous dimension matrix N <y AB is given by 

-Ln—l,n'—l ° J 



N AB 
JLn-l,n'-l 



(™-\) J Dwwr nl w^- l K AB { Wl ,w 2 ) . (6.47) 



in accordance with Eq. ( |6.40| ). Of course the same conclusions concerning the triangularity, 
universality and the eigenvalues of these anomalous dimension matrices can be drawn as in the 



10 The reader should be reminded that usually in the definition of the forward anomalous dimensions, Eq. ( 5.41 ), 
the value of n is shifted to n — 1. To choose the Mellin transform with the power (n— 1) rather than the power n 
is motivated by the interpretation of this quantity in terms of an angular momentum variable. However, choosing 
7i, as we have done here, leads to more compact representations. 
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non-singlet case above. However, the connection with the kernels K AB (w\,w%) by no means 
is trivial. However, computing the elements of the usual local anomalous dimension matrix 
j AB , individually, taking into account that only operators of the same value of N mix under 
renormalization, one obtains the (infinite) triangular matrices 



lZ> = J, dw * w * { 2 J, dw i wr*KFt?»u ^2)} , (6.48) 

7£ = 4n~-\y™. f 1 dw 2 w n 2 '- 1 {2 r W2 dw 1 w n 1 - n 'Kf(w 1 ,w 2 )\, (6.49) 



'0 I JO 

-1 , ( rl—W2 



= Kn'M«'/ dw 2 w^'l [ 2j o ~ W2 dw 1 wr n 'K^(w ll w 2 )Y (6.50) 
jf 1 dw 2 w%'- 1 {2 J^~ W2 dwi w^ n 'K^ G {w u w 2 )} , (6.51) 



GG _ (n-l\ (-) 

Inn' \ n i — \) u nn 1 



with 



(±) 



J(l + (-I)"-"') (l ± (-!)»'-***) = i(l±(-in(l± (-!)»'), (6.52) 



and 1 < n, n' < 00; a^, emerges in the pseudo-scalar case.0 The difference of the sign due 
to (— l) dB is compensated by an additional factor — (— l) dB which results from the different 
powers of w 2 in Eqs. ( |6.48| - |6T51"| ). The eigenvalues of these triangular matrices are given by 
their diagonal elements 7^f ■ In a way similar to Eq. ( |6.42 ) the splitting functions in the case of 
forward scattering are obtained for the singlet case [O, E9[ 



P AB (z)=2 dwO AB (w,z) / d Wl K AB { Wl ,w), (6.53) 



(no summation) with 



AB M = [ ]f W ~f I. (6.54) 

V 6(w — z)/z b(w — z)/z 



With respect to the cross terms an additional differentiation or integration has been introduced to 
remove the factors norl/n, followed by partial integration or changing the order of integration, 
respectively. The consideration of the pseudo-scalar case follows the same line of reasoning. 



6.3 Anomalous dimensions of twist— 2 operators at 0(a s ) 

After this general exposition we turn to the explicit calculation for the above operators 



Sj) in 0(a s ). We will consider first the unsymmetrized singlet anomalous dimensions K ( 



AB 



and AK AB for the unpolarized and the polarized case, respectively, in the a-representation. 
Here the parameters (0:1,0:2) are nothing but the usual Feynman parameters appearing in the 
calculation of the corresponding Feynman diagrams. As far as relations are concerned which are 
valid both for the unpolarized and polarized case we will, for brevity, only give the results for 
the unpolarized case in the following. The non-singlet anomalous dimensions obey in leading 
order 

K$ s = AK$ S = K q Q q = AK q q . (6.55) 



11 Note that according to our above convention on n the local anomalous dimensions are defined in the unpo- 
larized case at odd integers n and n', whereas they are defined at even integers in the polarized case in difference 
to other conventions in the literature. 
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For brevity we use the convention that the common factors resulting from the order of pertur- 
bation theory, as a s /(27r), and from the 0-structure of the different integration measures will 
be separated off from the anomalous dimension kernels. Then without introducing new symbols 
for the reduced kernels i^Q 9 (ai, a 2 ) we obtain 



K AB 



(i) {ai,a 2 ,K- 



2tt 



v «2 - «4, 



-d A 



(e(a 1 )9(a 2 )e(l-a 1 



a 2 )K AB ( ai ,a 2 ) (6.56) 



-l) dB Q(l - ai)0(l - a 2 )Q{a l + a 2 - 1)K AB (1 - a h 1 - a 2 



For the unsymmetrized kernels K A (ai, a 2 ) in leading order we get: 



1) Unpolarized anomalous dimensions 

if^ 9 (a 1 ,a 2 ) = C F |l - 8(ai) - 5(a 2 ) + S(ai) 



1 



a 2 



+ 8{a 2 ) 



+ §<J(ai)<J(a 2 )}, (6.57) 



i^0 G («l 5 «2) 

K^ 9 (ai,a 2 ) 
^ G (ai,a 2 ) 



- 2N f T R {1 — at — a 2 + 4«ia 2 } , 
-C F {8(a 1 )8(a 2 ) + 2}, 

Ca|4(1 — ot\ — a 2 ) + Y1ol\ol 2 + 8(a± 

1 



a 2 



2 + a 2 



(6.58) 
(6.59) 

(6.60) 



+ 8(a 2 ) 



2 + ai^|+ |/3 8{a 1 )8{a 2 ), 



where C F = (iV 6 2 - 1)/2N C = 4/3, T R = 1/2, C A = N c = 3, and the /^-function in leading order, 
Po = a — 4:T R Nf)/3. R The [ ] + -prescription is defined as 



/ dx[f(x,y)] + ip(x) 
Jo 



dxf(x,y) [ip(x) -tp 



(6.61) 



if the singularity of / is of the type ~ l/(x — y). 



2) Polarized anomalous dimensions 

AK™( ai ,a 2 ) 
AK^ (at,a 2 ) 
AK 3q (a 1 ,a 2 ) 
AK GG (a u a 2 ) 



K^(a x ,a 2 ), 

-2N f T R {l- ai -a 2 }, 
-C F {8( ai )8(a 2 )-2}, 

Y2C aol\Oi 2 - 



K GG [a x ,a 2 



(6.62) 
(6.63) 
(6.64) 
(6.65) 



Here (— l) ds has to be changed into — (— l) dB in Eq. ( |6.56| ). 

The anomalous dimensions for the polarized case have been derived for the first time in our 

[j]. Those for the unpolarized case have 
28f also the necessity of the 



13 and were later confirmed in Ref. 



previous paper 

been found several years before in Refs. |28|, [18] already. In Ref. 
symmetrization Eq. (|6.30|) has been pointed out. The kernels K AB and AK AB determine the 
respective evolutions of the operators 0™,O q 5 y and 0%, according to Eq. ( |6.28| ) in 0(a s ). 

As we have seen above the -^-representation reveals important properties of the kernels. It 
has been proven as a general representation in ||. The connection between the u>-variables and 

12 Note that the emergence of the term containing (3q is due to self-energy graphs only and does not induce 
scale breaking effects at this order, pPf. 
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the a-variables are given by Eqs. (|6.17|) . Using these relations we get the following expressions 
for the anomalous dimensions: 

Kf 1 f(w l ,w 2 ,K_) = ^^( K _)^-^(0(l- W2 - w Oe(l- W2 + ^ 1 )e(«; 2 )i? ( f i? («; 1 ,«; 2 ) 

+ (-l) dB &(l + w 2 -w 1 )&(l + w 2 + w 1 )e(-w 2 )K^ B (w u -w 2 )). (6.66) 
For the unsymmetrized anomalous dimension kernels in leading order we obtain: 
1) Unpolarized anomalous dimensions 

K q q (w u w 2 ) = C F \l-26(l-w 2 +w 1 )(l-- -) (6.67) 

I V [1 — W2 — Wi) + J 

-25(1 -w 2 - wi) (l - - -) ) + 3C F 6(1 - w 2 )5( Wl ), 

\ (1-W 2 + W 1 ) + J) 

Kf{w u w 2 ) = -2N f T R {w 2 + (1 - w 2 f - wj}, (6.68) 
K^ q (w u w 2 ) = -2C F {6(l-w 2 )6(w 1 ) + l}, (6-69) 

K GG (w u w 2 ) = C A \4w 2 + 3((l-w 2 ) 2 -wl) -25(l-w 2 + w l )(l + w 2 



(l-w 2 - Wx) + 

-26(1 - W2 - Wl ) (l + w 2 - - -) } + PqS(1 - w 2 )6(w 1 ). (6.70) 

\ (1 — W 2 + Wi )+/ ) 

2) Polarized anomalous dimensions 

AK qq (w u w 2 ) = K qq (w u w 2 ), (6.71) 

AK qG (w u w 2 ) = -2N f T R {w 2 }, (6.72) 

AK^ q (w l ,w 2 ) = -2C F {6(l-w 2 )6{ Wl )-l}, (6.73) 

AK^ G (w u w 2 ) = K GG (w 1 ,w 2 ) + ?,C A {wl-(l-w 2 f). (6.74) 

Again the factor (— l) ds in Eq. ( |6.66| ) has to be changed into — (— l) dB . 

Using Eqs. (CT-|g3T|, gg7|-g^), as well as Eqs. (gjjg-^jTgg, §TT]-g7§ in the pseudo- 



scalar case, we now derive the anomalous dimensions which are obtained in the local LCE directly 
from the non-local evolution kernels. Thereby the well-known relations 

1 1( ,n n-1 i 

I dw 2 - 2 — = E777 = ^ + 1 )+^ (6-75) 

(W 2 -1)+ 7^ni+l 



i a \b\ 

dw 2 w a 2 (l - w 2 ) b = ———— = B(a + l,b+l), (6.76) 
o (a + b + ly. 

are used. Here ip(x) = dlog(T(x))/dx, ^ F denotes the Euler-Mascheroni constant, and B(a,b) 
Euler's Beta-function. Again a common factor 

^ x2 ^ ^d g x 2a%, (6.77) 

in the unpolarized and polarized case will be suppressed, respectively. 
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1) Unpolarized anomalous dimensions: 



inn' 



inn 1 



Gq 
''Inn' 



7 GG 
I nn' 



1 ^ 1 

(„ + !)(„ + 2) +2 fy] 



(n + l)(n + 2) 



2 n' + l 
n — n' n + 1 



-C F 
C A 
+ 



n + l)(n + 2)(n + 3) 

' ~ 2 * 3n + A)5 nn , + 20 n 



(n 2 + 3n + 4) - (n - n') (n + 1) 



n(n + l)(n + 2) 
"1 2 
.6 



n+l 



n(n + 1] 
2n + 1 



= + 2V- + ^ 

(n + 2)(n + 3) ^ j 3C A 



j(6.78) 
(6.79) 
(6.80) 

(6.81) 



1 



,n(n + 1) 
with the following notation 



n — n 



(n-n' + 2) 



+ 



n(n+l) (n + 2)(n + 3) 



0, 



(±) 



|(l±(-l) n )(l±(-l)^ 

J 1 for n' < n, 
I otherwise . 



From this it becomes obvious that the symmetrization ( |6.26|) and ( |6.27| ) of the anomalous di- 
mension kernels only results into the property of the local anomalous dimension matrices being 
either odd-odd or even-even for the unpolarized or polarized case, respectively. 



2) Polarized local anomalous dimensions: 

' nn' 



11 
^inn' 1 



M qG ' 

Inn' 

M Gq ' 

Inn' 

M GG 

inn' 



-N f T 



n 



C 



(n + l)(n + 2)' 

[n + 3)8 r , 

4 



(n + l)(n + 2) 
1 
6 



-Or, 



n 

n+l 1 



+ 2V- + ^ 

(n+l)(n + 2) 3C A 



(6.82) 
(6.83) 

(6.84) 
(6.85) 



+ 



2n + 1 



1 



n — n 



7 ) - (n-ra' + 2)- 



n + l)(n + 2) 



As already noted in Section |6.1| the well-known forward anomalous dimensions are obtained as 
the diagonal elements of Eqs. ( p.78j - |6~8~T| , |6.82| - |6.85| ). 

The explicit computation of the anomalous dimensions is straightforward. In lowest order 
we have to determine the one-particle irreducible one-loop Feynman diagrams containing the 
considered operators as the first vertex. The new Feynman rules and an example of a calculation 
in the covariant gauge in the unpolarized case has been presented in |28| . Here we have performed 
the calculation in axial gauge which leads to essential simplifications in this order of perturbation 
theory. A sample calculation is given in Appendix O. 



13 In comparing these expressions with the values given in the literature 
[n, n') into (n — 1, n' — 1) sometimes has to be taken into account. 



33l the above mentioned shift of 



28 



7 Evolution equations and evolution kernels 



The distribution amplitudes and partition functions defined in Section 4 have to be determined 
experimentally at a given factorization scale /Zq. Their evolution in /i 2 , however, can be described 
within perturbative QCD and is ruled by renormalization group equations. In this Section we 
discuss generic evolution equations which directly follow from the renormalization group equation 



7.1 General properties of twist— 2 evolution equations 

Let us first discuss general properties of the evolution equations and their kernels at any order of 
perturbation theory. For simplicity we restrict the discussion to the non-singlet (NS) case since 
the arguments in the singlet case are quite analogous. 

The evolution equations of the partition functions in question are consequences of the renor- 
malization group equations of the complete twist-2 operators, i.e. including the trace terms which 
are proportional to x a . These operators are summed up local operators of definite twist. Despite 
the fact that the vector operator 0* wlst2 and the scalar operator O twlst2 = x (r O^" st2 (Kix, K 2 x) 
fullfil the same renormalization group equation (|6.1|) , as has been shown above, the evolution 
equations of their matrix elements are different. 

For general values of the arguments Ki, k 2 the operator matrix elements read 

(P2|0 W2 (kiz,k 2 x)|pi> = ifj(K_x P+ ,K„xp^,K 2 _x 2 ,p m ,fi 2 )e lK + xp -Q J (x-p 2 , Pl ) (7.1) 

and 

(p 2 \0^' twist2 (K 1 x, K 2 x)\pi) = iQ Jll (p 2 ,p 1 )Fj(K_xp + , K + xp., k 2 _x 2 , Pl p 2 , /i 2 ) (7.2) 

+ iQ J (x;p 2 ,p 1 )i {k^8 k _ xp _ + K + p^_d K+xp _ + K_p%d K+xp _ + 2x fl K 2 _d K 2_ x2 } 

XF/, (/€_££>+, K^XP-,K 2 _X 2 ,Pip 2) [l 2 ) 

which is nothing but the Eq. ( |4.1(J| ) in the coordinate space. 

Whereas for the representation of the Compton amplitude the operators and their matrix 
elements can be taken at arbitrary values of x, cf. Section 5, for the renormalization properties 
considered in Section 6, it is essential that the operators are defined on the light-ray x, x 2 = 0. 
Furthermore, in the subsequent considerations the operators have to be taken on the light-ray 
at general values of Ki . 

The matrix elements of the scalar operators N obey the renormalization group equation 



2 d 



dfr 



ip 2 \0 JiS {K 1 x,K 2 x)\p 1 ) = I d 2 K/-f NS (K 1 ,K 2 ;K' 1 ,K' 2 )(p 2 \O m (K' 1 i,K f 2 S:)\p 1 ) (7.3) 



and correspondingly for 0^ s with / ~f^ rs = A7 7VS . Having in mind the kinematic decomposition 
of the matrix elements, Eqs. ( |4.4j , |7.1|) we get renormalization group equations for the partition 
functions in coordinate space, fj, separately: 

V 2 ^fj(K-xp +1 K_xp^y K +* p - = [ d 2 «V 5 («i,«2;«i,«2)/j(«'-5p+>«'-5p-)e te ' + *- • (7.4) 



dy 2 

In the same way we obtain the renormalization group equations for the matrix elements of the 
twist-2 vector operator ( |7.2[) and the corresponding partition functions Fj. After the kinematic 



29 



decomposition the operator relation Eq. (3.8) implies 



d 
dfi' 



;Fj(K_Xp + , K_Xp_, K + Xp_) 



c? 2 k'7 (ki, k 2 ; k[, k' 2 ) Fj(k'-XP + , k'.xP-) . 



By differentiation w.r.t. xp_ and xp + the equations 



// 



d 
djj? 



_d_ 

djj? 



K-d K _ x 



xp + 



Fj 
Fj 



dV 7 k 2 ; k[, k' 2 ) n-d K _ xp _ + K+d, 



K+Xp- 



d !Sl 7 ( K i) K 2! K 'n Wjj) k '-^k' 



(7.5) 

^ (7.6) 
(7.7) 



are obtained with the same arguments as in Eq. (|7.5|) . The renormalization group equations are 



valid in the limit x 



x, i.e. x 



0. Therefore a further relation 



d 
djj 2 



^-d K 2 x 2F(n_xp + , K_xp_, n + xp^, k 2 _x 2 " 



(7.8) 



d 2 K j ns (ki, k 2 ; nf t , k' 2 ) K /z _d K ,2 x 2F(k' _xp + , n'_xp_, K,' + xp_, k' 2 _x 2 ) 



holds. Here we will not discuss this special evolution equation in detail. We remark however 
that the corresponding anomalous dimension belongs to the operator 



0'{k\x, K 2 x) 



X 



d_i, 
dx~V 



(/«ix)(x7)C/(ki, k 2 )iP(k 2 x) - iP(k 2 x)(x^)U(ki, K^ip^x) 



(7.9) 



and has to be calculated independently. Here x* denotes a second independent light-like vector 
which obeys x*x = 1. 

Note that there are two evolution equations, (|7.4j) and ([7.5|) , for different partition functions 
which are non-trivially related by 



Fj{K_Xp + , K_Xp_, K + Xp_ j 



d\ fj(XK_xp +1 XK_xpJ)e 



iXn+xp- 



(7.10) 



We finally show that this relation is compatible with applying the renormalization group oper- 
ator. The evolution equation for Fj reads 



2 d 



dXfj{X^xp + ,X^xp^)e lXK+ip - = /dV / dA 7 NS (A/t, rt)fj{rt k _xp^)e iK 



djj? Jo 

in the representation of Eq. ( |7.10|) . Changing the integration variables k' = Xk and applying the 
scaling relation Eq. flb.ll ) for the anomalous dimensions 

7(Ak,Ak) = X- 2 ~f(K,k) , (7.11) 



we obtain 



/'" 



_d_ 

djj? Jo 



dX fj(Xn_xp + , A/t_xp_, Xn+xp. 



d 2 «7 N (k, Bl) I dX fj(Xk_xp + , Xk^xp-)e 



iXh+xp- 



which is the evolution equation (|7.5|). This shows that Eq. ( |7.5| ) and Eq. (|7^|) are equivalent. 
The evolution equations ( [7.4| - [77^ ) yield already a complete description in the non-singlet case. 
However, commonly the evolution of parton densities or distribution amplitudes is described 
in terms of functions which depend on the related momentum fractions. This description is 
obtained after a Fourier transformation of the above functions to the functions f(z + ,Z-) and 
F(z + , Z-) which have been introduced in Section 4 before. In the following we will consider the 
evolution equations of these quantities. 
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7.2 Two- Variable Equations 

In this paragraph we derive the evolution equations for the double distribution amplitude in the 
singlet case. Let us start our consideration with the double-variable distributions fj\z+, zJ) 
(z + , Z-). In accordance with Eqs. QOj) and (|4.22| ) we use the kinematic decompositions: 

K-fa\0*^ 2 (-K-X,K-5;\px) = (xp + )- 1 ff ) (K^p + ,K^,p lP2 , l x 2 )Q J (x;p 2 , Pl ), (7.12) 

K 2 _(p 2 \0 G ' twiSt2 (-K_X,K_x)\ Pl ) = (xp + )- 1 f^ 2 \K^p + ,K_Xp.,p 1 p 2 ^ 2 )Q J (x ] p 2 ,p 1 ). (7.13) 

Combining both cases we use the notation f A , where A = q, G, thereby suppressing the index 
(2). We are interested in the evolution equations for the physically relevant Fourier transforms 

#(*+,*-) = r K-iW-)e^ S ^ J »-'- ) , (7-14) 

J-oo Z7T J-oo Z7T 

where variables being irrelevant for the present considerations have been omitted. In order to get 
the evolution equations for the singlet case in z-space we have to perform a Fourier transform 
of the equations analogous to (7!3). We take matrix elements of the renormalization group 
equations (|6.29 ) with the corresponding evolution kernels 

K AB ( Wl ,w 2 ,K-) = K d _ B - dA K AB (wi,w 2 ) (7.15) 

and perform the Fourier transform to get the evolution equations for f A . Inserting the definition 
(|4.14j ) into the right hand side we obtain 

f—jr\^-) = r dK ^ p+ r dK -w- e w + -^- e <K-^+*-*-) (7.i6) 



x J DwK AB { Wl ,w 2 ,K_)K d _ A J Dz' (^y d B e -i-'-HP^' + +P-^)fBm^ 
Carrying out the integration over d(K-xp±) we get 

» 2 4- 2 ff [2 \^Z-) = jDwK AB ( Wl ,w 2 )w 2 d - 



x[dz'_dz' + 5(z^-w 2 zi+w 1 )5(z + -w 2 z' + )ff (2) (z' + ,z'_). (7.17) 



dfi 2 



This finally leads to 

» 2 -^- 2 ft (2) (z + ,Z-) = jDz'T AB (z + ,z.,z' + ,z'_)f? i2 \z' + ,z'_), (7.18) 

r AB (z + ,z_,z> + ,z'_) = J ^-K AB (w 1 = -z_ + Z -±zL,w 2 = Z -±) ( Z ^) dB . (7.19) 

\z + \ \ z' + z + J \z+J 

Now we consider the kinematic decomposition for the first parametrization: 

(p 2 \0^ twist2 (-K„x,K„x)\ Pl ) = iff\K„xp + ,K^p_,p lP2 ^ 2 ) Q J (x;p 2 , Pl ) (7.20) 

(p 2 \0 G ' twmt2 (-K_x, K-x)\pi) = i(ixp + )f ( f {1) (K_xp +1 K_xp^,p 1 p 2 ,fi 2 )Q J (x-,p 2 ,p 1 ). (7.21) 

The derivation of the evolution equations for the distribution amplitudes Eqs. ( |7. 20|j7T2l|) follows 
the same line. The only difference is that the integration with respect to does not lead 

directly to the 5-function and one obtains 

» 2 /-J?(z + ,Z-) = J Dz'T AB W(z + ,z.;z' + ,z'_)f B (z' + ,z'_), (7.22) 



31 



with 

T AB ^\z + ,z.;z' + ,z'_) = J DwO AB (z + -w 2 z' + )K AB (w 1 ,w 2 y(z^-w 2 z'_+w 1 ), (7.23) 

(no summation), where AB is given by 

, f 8(z) for cIb — cIa = 

AB (z) = — / e ^-a P+ , = J 9 S ( Z ) for d _ d = J (7 24) 

Here denotes the sign-function, cf. Eqs. (|8.15| , |D.4j) . 



2 e(^) f° r ds — d,A = — 1- 



7.3 One- Variable Equations 

The equations given above cover the most general case. These equations depend on two partition 
variables z + and Z- . The question arises whether it is possible to express them in terms of single- 
variable distributions imposing a one-parameter constraint. Having (physical) scale invariance 
in mind it is reasonable to introduce a kinematic condition by 

r = ^ , (7.25) 

xp + 



under which the evolution equations are derived. This approach has been outlined in H 13 . It 
is possible to define the matrix elements formally by 



(p 2 \O q (-K_X, K_x)|pi) 



(ixp+) 

(p 2 \0 G (-K^X, K-X)\px) 



(ixp^ 



xp—=Txpj r 



/+oo 

dte-^-^+^q^r), (7.26) 
-oo 

[ +OC ' dte- lK - ip+t t$ G {t,r) . (7.27) 



Here $(£, r) denotes, in the above sense, a partition function which in the limit r — > — 1 tends to 
the standard Brodsky-Lepage wave function. Note that the parameter r is related to x-space. 
This parameter may change into the purely kinematic variable r\ = qp_/qp + in some applications 
after a Fourier transform. In Eq. ( p.29| ) 

/ Dz5(z + - t + V z-)ff {1) (z + =t-r,z., z.) = ff l \t, V )t dA ~\ (7.28) 
7] occurs naturally as a variable in momentum space. The reverse formula reads 

, _j_ / +o °* r 00 ^ i-i^r,-(^-') u-iw, 



^"(— ^^exp^-^f^-Vf"^,,). (7.29) 



Eq. (|7.29 ) has been derived by assuming, that the restrictions of the support are completely 
included into the definition of the functions fj (t,rj). This can be achieved using the repre- 



sentation ( |3.15| ) for the Heaviside-functions and performing the limits lim e ^ + later on, cf. [53 
One may equally well consider the set of one-variable functions ff (t, rf) instead of the two- 



variable functions fj\z + ,z-). Although the Compton amplitude depends on the functions 
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Ff^\z_, z + ) our evolution equations for simplicity have been formulated with the help of the 
functions ff^\z-, z + ). However, both are related by 

F?% + , *-) = [ § ft m (i, y) 6(A - |, + |)6(A - |z_|). 
The relation for ff corresponding to Eq. ( |7.28|) reads 

Ff^it,^- 1 = J DzFf 1 \z + ,z_)6(z + -t + r ] Z-), 

which yields 



(7.30) 



1 d\ 



■A{1) 



t 



d A -l 



(7.31) 



We finally obtain the one-parameter evolution equations for the functions fj and / 

d — r+oo 



Mi). 



,,2 " \ 

djj? {t ' T) 



-oo 
+oo 



dt'V^(t,t f ,r)^ s (t',r), 
dt'V^(t,t',r)f^\t',r) 



(7.32) 
(7.33) 



with fj S — ff — ff . The corresponding extended kernels V^(t, t', r) read 



dati 
o Jo 



l—ai 



da 2 K (ax,a 2 ) 



d(nxp-i 
2V" 



d B -d A 



— exp [inxp + t — (1 — ol% — a 2 )t' + r(«i — a 2 ) } . (7.34) 



They obey the scaling relation 

V£M,t) = V^{t,t'-r) = -K4 S (-,-,l). (7.35) 

r \t T / 

For convenience we rewrite the general expressions for the evolution kernels in the variables 
In one-loop approximation they are given by 

Kxt (*, t>V= n _ V t {t, t , T) 



2tt 



(7.37) 



where 



V q Ut,t' } r) = l-{v^(x,y)p(x,y) + V^(x,y)p(x,y) + IC F 5(x-y)} (7.38) 



V %(t,t',r) = — 



1 /2y-l 



2r V 2 
1/2 



{v^ G (x,y)p(a;,y) - ^ G (x, y)p(x, y)} 
V %(t,t',r) = ±(-^^{v G \x,y)p(x,y)-V G \x,y)p(x,y)} 



(7.39) 
(7.40) 
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K£t(M',r) 



2y 



2r \2x 
Zt I 



^ {V GG (x,y)p(x,y) + V GG (x,y)p(x,y)} (7.41) 



A*3L(t,f,r) 
AV^t^r) 

AV^(t,t',r) 

with the O-structure 



1 /2y-l 



2t 
1 



2r V2x - 1 
1 /2y-l 



2r V2a; - 1 
+ 7TnPoS(x-y) 

ZT I 



{AV qG (x, y)p(x, y) - AV qG (x, y)p(x, y) } 
{AV Gq {x, y)p{x, y) - AV Gq (x, y)p(x, y) } 
{ AV GG (x, y)p(x, y) + AV GG (x, y)p(x, y)} 



(7.42) 
(7.43) 

(7.44) 

(7.45) 



p(x,y) = 61- 



x 



e 



X 



sign(y) , 



(7.46) 



The non— singlet kernels are given by Vg^ t — Votxt- 
different partial expressions V AB and AV AB read 



and AVj 1 ^ = AV^ q ext , respectively. The 



V qq {x,y 
V qG (x,y 
V Gq (x,y 
V GG (x,y 

AV qq (x,y 
AV qG {x,y 

AV Gq {x,y 

x,y 



AV GG 



2N f T R - 

y 



l l 

y (y-x). 

X " 



4(1 -x) + 



1 - 2x 



x 



2-, 



C F 1 - 
L 1/ 



2/ 



-2N f T R - 

y 2 



3-2x + 



x 



1 — X 



+ 



(y - x). 



c F 

C A 



If 

x 2 

y 

x 2 
2 - + 

y 



y + x 



(y-x). 



(7.47) 
(7.48) 
(7.49) 
(7.50) 

(7.51) 
(7.52) 

(7.53) 

(7.54) 



Details on the [ ] + -prescription, which acts to the right, are given in Appendix [C|. Note that 
the kernels given in Eqs. ( |7.38| - |7T45|) apply to the full range of variables, i.e. they represent the 
kernels completely. The function V qq t (t,t' , r) was already derived in Refs. [35], Q. 

The determination of these kernels is straightforward observing, however, that these quan- 
tities are distribution-valued. In the calculation one may set r = 1 and later apply the scaling 
relation, Eq. (|5.11| ), taking into account the general 0-structure in the (t, t')-plane derived in 
0. Alternatively, the calculation can be done in the different kinematic regions separately. A 
sample-calculation for the kernel V Gq is given in Appendix |D] All other known representations 
in the literature consist explicitly of two or three parts being calculated separately. 
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cf. 



Recently, the kernels (^3|-|^4l|) were also calculated in Ref. ||11|| , using a different notation, 
for related work. For r = 1, corresponding to £ = 1 in they agree with the results 
A further independent representation has been given recently in Ref. 



obtained above 

parameter r introduced above equals to — £/2 in the notation of |37| 



37[. The 



8 Special cases 

The evolution kernels given above cover a series of limiting cases which were studied in the 
literature before. These are characterized by special kinematic constraints for the processes, as 
in the case of forward scattering (p 2 \ — > (pi\ = (p\ or the transition from the vacuum state (0| 
to a hadron state {p\ being related to the hadron wave functions. 



8.1 The Brodsky— Lepage limit 



For r = ±1 the equations ( |7.38| |7T45]) transform into the limit (p2\ — > (p\, (pi\ — > (0|, which 
is known as the Brodsky-Lepage [38] and Efremov-Radyushkin [|39| case. This limit may be 
performed formally leaving p\ — > 0, which leads to correct results, cf. j3"5| . The corresponding 
evolution equations are obtained using as variables x and y, Eqs. ( [7.361 , |7.32| , |7.33| ). As an 
example we consider the simplest case. For < x, y < 1, |r| = 1 one obtains 



V*"(tJ) = \c F {Q(y-x) 



x 1 

y y (y 



l 



+ Q(x-y) 



I — x 

^y 



l 



l 



y + (x- y). 



+t;6(x 



/)}■ 



5.1^ 



8.2 The 'near— forward' representation 



This representation has been introduced by X. Ji and has also been given in Ref. ||11|| . It contains 
the forward case as limiting case. A correct application of the evolution equation for near forward 
matrix elements needs the representation for the region t < r and t > r separately. They are 
obtained from our general kernels as special cases. As an example we show here that our general 
structure for V^t covers also the case t > r, i! > r. Note that 



sign y = —sign y and 0(1 — x/y) = Q(y — x) 



(8.2) 



One therefore obtains 



V qq ext (t,t',r) = \c F Q{y-x)\- 
= \c F Q(y - x) 



1 + 



(3/-2O+ 



(y - x). 



1 + 



yy 



1 — X 

^y 



1 + 



ix-y). 



} + p F 8(x-y) 



+ jC F 8{x - y) 



(8.3) 



for r = 1. This representation was also given in Ref. JTT 



14 Several independent calculations of the evolution kernels for the meson wave functions were performed in 
Refs. pOj. 
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Starting from the anomalous dimension we reproduce the results of |37] and generalize them 
even to the region with t' ^ 1. We obtain 



K GG (t,t>,0 



AK qq (t,t' } £) 

AK Gq (x,0 
AK GG (x,0 



C 



t 2 + f 2 - e/2 
F it' 2 -eim'-t) + ' 2 



C, 



(f 2 - e/4) 2 

t' 2 + (f - t) 2 - e 2 /4 



t 7 



2C 



t (t> 2 - e/iy 



if 2 - 
(f-t) + 



-t(3t /2 -e 2 /4)-(t' + t)(t'-t) 5 

^ 9 (M',£), 

t 2 - (t - tf - £ 2 /4 . 



+ t'(t' 2 + £ 2 /4) 



it' 2 - e/4) 2 
t' - (f - t) 2 - £ 2 /4 



t(i' 2 - £ 2 /4) 



2 % (f 2 -e 2 /4) 2 



it' 2 - £ 2 /4) 5 



+t'(t' 2 +e/i). 



-t(3t' 2 - £ 2 /4) - 2t'(t' - tf 



+ ^S(t> - 1) 



(8.4) 
(8.5) 
(8.6) 
(8.7) 



(8.8) 
(8.9) 

(8.10) 

(8.11) 



To get a complete representation of the kernels the independently calculated parts for t < — £/2 
and \t\ < £/2 have to be added, which can be directly derived from the general expressions given 
above. 



8.3 The Altarelli-Parisi limit 

Let us consider the case of forward scattering p 2 = p\ = p. The corresponding evolution kernels 
can be obtained after some calculation from Eqs. ( [7.38| - [7T45] ) in the limit r —>■ 0, if > 0, 

KmV f xt (t,t',T) = ±P AB (z)Q(z)Q(l-z) with z = t, (8.12) 

cf. Appendix [D]. Alternatively, the splitting functions for deeply inelastic forward scattering can 
be obtained by direct integration of the kernels Eqs. ( |6.57| - |6l)5|) , which is performed in the 
following. The splitting functions are obtained by 

P AB {z) = P AB (z)6(z)6(l-z) (8.13) 
du9(l-u)9(u)(l-u)0 AB (u,z) / d£^ B ((l-^)£, (!-«)(!-£)), 

-oo JO 



where 



6 AB (u,z) = ( f~ U \, d R f Z -$ )• (8.14) 
1 —6{z — u)/z o(z — u)/z 



15 In ref. ||37| t he variab le t was denoted by x, having a different meaning in our notation. Therefore we use t 
instead in Eqs. Here our variable r has to be identified with the variable — £/2 from |37|]. 
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The 0-structure resulting from the measure Da is universal. In deriving Eq. (|8.13 ) it is useful 
to apply the relations 



9{x) 



lim — 

£^0+ 27T 



ii + e' 



1 

2^ 



(life 



(8.15) 



which are valid for tempered distributions [p4| |. P AB (z) (AP AB (z)) are the well-known splitting 
functions [33| for unpolarized and polarized deep-inelastic forward scattering in 0(a s ): 



P qq (z) 
P qG (z) 
P Gq (z) 
P GG (z) 

AP qq (z) 
AP qG {z) 

AP Gq {z) 
AP GG (z) 



1 + z 2 



2N f T R 



1-z 

2 



z 2 + (l- 



^2 



C, 



1 + (1 



Z 



2C A 

P gq (z), 
2N f T R 

2C A 



z 

1 1 



2 + z(l-z) 



2z + 



(1-*). 



(8.16) 
(8.17) 
(8.18) 
(8.19) 

(8.20) 
(8.21) 

(8.22) 
(8.23) 



The parton densities of the quarks q(z) and antiquarks q(z) are given in Eqs. ( |5.22| , [Slffi ) above, 
with z denoting the momentum fraction carried by the partons. G(z) denotes the gluon distri- 
bution. Their evolution equations read in leading order for the unpolarized case 



d 



g NS (^ 2 ) 



a s (/i 2 ) f 1 dz' 



dy? 
(I ( S(z,/i 2 



2tt 



z' 



dn 2 V G(z } /j, 2 



a s (fi 2 ) f 1 dz'^~ ( z 



2tt 



— P - 



7K b c*V), 

G(z',/i 2 ) 



(8.24) 
(8.25) 



Here g NS (z ;/ u 2 ) denotes a non-singlet combination of quark densities, as e.g. q(z,fi 2 ) — q(z,fi 2 ), 



T,(z,fi 2 ) = X^2i qj(z,fi 2 ) +qj(z,fi 2 ) is the flavor-singlet combination of the quark densities, 

and P is the matrix of the leading order singlet splitting functions. The respective evolution 
equations in the polarized case are obtained by replacing q, q and G by Aq, Aq, and AG and 
the splitting functions P AB (z) by AP AB (z). 



9 Solutions 



Here we extend the method used in Ref. JT0| for the non-singlet case to the singlet case. The 
idea follows the original solution of the Brodsky-Lepage equation given in p9[. We start from 



Eq. ( 7.17 ) and form the moments of the distribution functions 



-n 



dz + zlf A M(z + ,z_) . 



(9.1) 
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leading to 



with 



(9.2) 



r£ B ( z _,z'_) = J Dw5(z^-w 2 z , _+w 1 )w^ dB K^ B (w u w 2 ) 

= \f Q dw 2 B(l + z_- w 2 (l + z'_))Q(l -z-- w 2 (l - z'_)) 

x0(l -z'_)Q(l + z'_ )K$ B {-z_ +w 2 z'_,w 2 )w^ dB . (9.3) 
In this way a first diagonalization is obtained. The kernels r^ B (z_, z'_) read in the unpolarized 



case 



+ Q(z' - z) 



1 + z'J \n (z' — z) 



35(z-z'), 



(9.4) 



-2N f T R U(z - z>) ^ + 2 < Z - - ^' - * 



+ Q(z' - z) 



(n + l)n(n — 1) 
l + z\"n 2 + 2n(z' - z) - (2zz' - 1) 



(9.5) 



^ q (z,z') 

r° G (z,z>) 



.1 + z'J (n + l)n(n - 1) 

c A {e {z - ~J) (l^X (W + M-n + ^ + 3(i + ^)d-^) 

[ VI — zv \ n— 1 ra n+1 

2 2 \ 

+ 



n (z — z') + J 

+ Q( z > - Z ) (i±£)" p(l-^)(l + ^) + 6^ + 3(1 - *')(! + *) 



2 

— + 



n- 1 
2 



Similarly one obtains for the kernels in the polarized case 



n n + 1 

+ p S(z-z'). 



(9.7) 



ATf(z : z f ) 



C F \Q(z-z') 

+ e(z' - z) 



l-z\ n (l 2 
■ ' + 



1 — z' ) \n (z — z'). 
l + z\ n fl 2 



1 + z'J \n (z' — z). 
1 -z\ n 1 



+ 35(^ - 



(9.8) 
(9.9) 



c F {e {z -z>)(— ) - + e( ,<-^— j J, (9.io) 



l-^\ n 1 „, , s /1 + Z\ n 1 



atJW) = -^ /ra {e(,-y)(^)"i + e ( ^-,)(i±l)"i} 
Ar^(,,,') 



+ 6(/ - ^) 



1 - z 
l-z' 
( 1 + z\ n (4 



4 2 
+ 



n \z' — z) 



+ P S(z-z'). 



(9.11) 
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As a next step we have to find a solution of Eq. ( |9.3|) . One finds that the following symmetry- 
relations hold for kernels T GG (z, z') and AT GG (z, z') : 



y /2\n-pAB 
- ) 1 n 
Jl\n a WlB 



Ar* B (z',z)(l-z 2 r ■ 



(9.12) 
(9.13) 



This is expected because of similar problems studied previously in |10|, |39[ . The final diagonal- 
ization of Eq. (|9.2| ), referring to the partition functions f A ^ 2 >(zJ) for which the quarkonic and 
gluonic operators are dealt with equally, can thus be performed using Gegenbauer polynomials 



Cl(z), cf. [41 



AT AB 



AB 'z,z') 



ECi 1 - ^Cl^z)C^{z')N-) n 



z, z 



E Ar^(i - zrc n m +1 '\z)c^i\ z ')N-l . 

The coefficients N n>m are the normalization factors of the Gegenbauer polynomials 

r+l 



(9.14) 
(9.15) 



dz(l - z 2 TC? +1/2 (z)C n k +1/2 (z) = S lk N n , k 



6. 



ik 



2 2n + i r 2/ + 1\ (k + n + 1/2) ■ k\ 



i -i 



2 J 2irT(2n + l + k) 



(9.16) 

T A ^ and AT A ^ are the respective expansion coefficients which can be easily calculated and are 
found to be the diagonal elements of the triangular matrices discussed in Section 6. 

Another method of solution consists in forming also Mellin moments in the variable z_, 



-i 



dz.z k _f^\z.) = f$> . 



(9.17) 



The evolution equations may be written as, cf. also |28 

o d 



// 



d^ Jnk 

with the transformed kernels given by 



,A(2) _ a S {^ 2 ) y> V AB fA(2) 
J nk r>_ / i n.kl.J nl ) 



2tt 



(9.11 



1=0 



n,kl 



Dw + w^ +l - dB w k 1 - l K AB {w 1 ,w 2 ) 



k 
I 



(9.19) 



These evolution equations are not diagonal with respect to the indices (k, I). The explicit ex- 
pressions for T AB t and AT A kl are given in Eqs. ( |6.7S| - |6~81"| , |6.82| - |6.85|) . For fixed n they form 
triangular matrices. The eigenvalues are the diagonal elements k = I 



r& = / Dw + w n 2 +k ~ dB K AB ( Wl , W2 ) = ^ B k _ dE 



(9.20) 



The coefficients j AB k _ h . are the anomalous dimensions of the forward case with a shifted Mellin 
index. A solution for the singlet evolution equations in the case of the single-variable equations 



was given in [42 recently. 
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10 Conclusions 



The scaling violations of the various observables, which occur in space-like processes, can be 
traced back to the renormalization of light-ray quark and gluon operators, which define these 
quantities. Therefore one may obtain a general description by calculating first the anomalous 
dimensions, or evolution kernels, of these operators and establish their connection to the different 
observables through the respective integral mappings and, in some cases, subsequent limiting 
procedures, afterwards. In the present paper we were limiting the investigation to the case of 
twist-2 operators and massless QCD. They are found by a corresponding twist-decomposition 
starting with a genuine operator representation. The space-like evolution kernels are process- 
independent. To form the observables one has to calculate also the process-dependent coefficient 
functions, which are trivial in lowest order. Furthermore, expectation values of the operators 
have to be formed, allowing for the description of a wide variety of processes. 

In the light-cone expansion both vector operators O a and scalar operators x a O u appear. 
As was shown they are equivalent under renormalization and it suffices therefore to study the 
evolution of the simpler scalar operators. The respective relations for the vector operators can be 
obtained by integral relations. This holds also for the relation between their general non-forward 
expectation values between the states (pa| and This connection is an important one, as it 
forms the basis for the integral relations by Wandzura and Wilczek [[HI and one of the authors 
and Kochelev ||20|| . These relations were derived earlier in the local operator product expansion 



and a subsequent analytic continuation from the integer moments and a Mellin-transform. Here 
it is obtained directly by the integral relation between the scalar and vector operators. 

One may derive a series of all-order relations for the evolution kernels of the light-ray opera- 
tors. In the local representation they are given by infinite triangular matrices 7^ B n2 , the diagonal 
elements of which are the anomalous dimensions in the forward scattering case. The anomalous 
dimensions are independent of the the total spin n\ + HI2. 

The expectation values of the operators between the states (p2 1 and \p\) are complex-valued 
distribution amplitudes through which the scattering cross section may be represented, which 
are, however, no observables. Nonetheless they obey evolution equations which are implied by 
the renormalization group equations of the operators. Commonly these matrix elements are 
represented by a two-fold Fourier transform as distribution amplitudes which depend on the 
two momentum fractions z + and z_ along the momenta p + and p~. This representation is a 
generalization of the case of forward scattering in which only one momentum fraction emerges. 

The two-variable amplitudes may be related to single-variable amplitudes by implying a 
kinematic constraint to the general case, i.e. demanding r = xp_/xp + = const. In this represen- 
tation the kernels depend on the variables t and t' besides the parameter r. We obtained 
a unique representation for the whole kinematic range of t and t' . One may easily describe the 



scaling violations in the case of vacuum-meson transitions [[38], |39| in the limit r — > — 1. With the 



help of the functions V^f(t, t', r) other recent representations [|Tl|, |37j can be verified after some 
calculation. The evolution kernels in the limit of forward scattering, which are the well-known 
splitting functions, can be obtained by an integral relation form the two-variable kernels at one 
hand. They result also from the single-variable representation in the limit r — > for t' > from 
the kernels V e ^(t, t', r). 

The solution of the two-variable non-singlet and singlet evolution equations require a diag- 
onalization of the evolution kernels with respect to the variables z + , z' + and z_, z'_, respectively. 
Here the w-representation shows, that the diagonalization with respect to the first variables can 
be performed by a Mellin transform, as in the forward case. This is an all-order property. The 
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diagonalization with respect to the second set of variables through Gegenbauer polynomials re- 
lies on a symmetry property of the kernels, which exists at leading order, but is broken in higher 
orders. Alternatively, one may perform a second Mellin transform. The latter representation 
is, however, only suited for discrete representations over integer moments. An easy analytic 
continuation is lacking in the latter case since the eigenvectors do not obey an analytic relation 
which can be continued to complex numbers. 
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A Kinematic Relations 



The Compton process is characterized by the invariants 



(Pi + <?i) 2 
(Pa - Pi) 2 



(P2 + 92) 
(92 - <?i) 



g? + 2 Piqi -- 

Ml - 2 Pl p 2 



Ml 



Ql 



2^292 



-Ql -Ql- 2 gi q 2 



(A.l) 



Here M 12 are the masses of the initial- and final state hadron. q\ 2 — ~Qi 2 are the virtualities 
of the incoming and outgoing photon. The initial-state photon is assumed to be space-like 



Ql = -ql > 



(A.2) 



while the final state photon can either be light- or time-like with 



2 „ 2 
Q2 < ™*, 



(A.3) 



to avoid s-channel mass thresholds in the subsequent description, as they occur in reactions like 
7* + p — > V + p', where p and p 1 are the initial and final-state hadrons, respectively, and V 
denotes a meson at the photon-side. In the latter case the process is no longer space-like and 
the light cone expansion cannot be applied. On the other hand, one may produce mesons on the 
hadronic side in a space-like process. 
The inelasticity v is given by 



p +qi + -{Ml - M; 



and the cms-energy s reads 



p + q 2 -^{Ml-Ml) , 



(A.4) 



2z/+ - 
2 



ql + Ml + M$-t 



(A.5) 



The virtuality Q 2 



-q 2 is 



Q 2 



t 
4 



\ {Ql + Ql) 



(A.6) 



We define 



The variable £ is given by 



x 



Ql 



V 



(A.7) 



x t + 2Q 
4 ~ 2 4z/ 



4^ 



(A. 



showing the relation of the variables £ and 77. For tjv — > 0, as assumed in the generalized Bjorken 
region, their difference is x. In the case of forward scattering, rj = t = 0, £ and x turn into the 
Bjorken variable xb = Ql/i^P-Qi)- Also other kinematic domains were studied for the Compton 
amplitude, as e.g. wide-angel scattering p\.p 2 — > 00 j43|, a quantity being fixed in the present 
treatment. 
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B Sample calculation of the anomalous dimensions 

In this Appendix a sample calculation for the non-local anomalous dimensions is performed. As 
an example we will consider the derivation of the anomalous dimension AK Gq . The calculation 
is performed in the MS-scheme using dimensional regularization in N = 4 — 2e dimensions. The 
anomalous dimensions are given by the pole parts (p.p.) 

I=\ /p.p. + 0(e°) (B.l) 
of the respective integrals. The gluon operator, Eq. ( |6.8| ), reads in the momentum representation 



Of(«i,«2) = \x^' F a ^{K l x)F%, v {K 2 x)-F a llv {K 2 x)F a /{K l 



X 



dpidp- 



4(2tt) 



2 / IKlXpi— lK2Xpi— IK1XP2 ~> 

2N ^ C ' 



V,^P 1 ^(Pi)(p^(P2) -P^(p 2 )) . (B.2) 
The one-loop correction to this operator is calculated. We work in the axial gauge, so that 

X^Fp V \KX) > xd K x-A. v yKX^ 

x^F^^kx) > x^£^ va pd K ^.A p ( K Kx) . (B-3) 

We are interested in pole-parts only, which for AK Gq correspond to the counter term of the 
quark operator. Actually not the complete gluon operator but only the term 

x^'f;( Ki x)F^(k 2 x) (B.4) 

needs to be used in the calculation which corresponds to the analytic expression: 

d N Pl d N p 2 1 —r. 



T^X^ F^{k 2 X)FJ{k X X)S) = J T^yv ^bl)^^!, «2, Kl)^(P2) , (B.5) 



with 



and 



/(p 2 ,Pi, «2, «i) = e^-^^JiCpi,^, «!, «a) + e*™-*™/ 2 Cp 2 , Pl , k 2 , «*) (B.6) 



/" k 1 ( 

h = J d N k(igj tl )i-^(ig'y v )-r lg 



lin _ x^ + x a q^ \ ]_lf up _ x»q p 2 +xPq» 2 \ 1_ 
xqi ) q\ % \ xq 2 J q\ 



x (-iq1x x i{xq 2 )e KXpa ) e lKX \ (B.7) 

where q\ — pi — k and q 2 = p 2 — k. I 2 (p 2 ,pi, k 2 , k{) is obtained by interchanging •y u <-> 7^, 
Pi ^ P2 and «i <->• k 2 . Here contrary to our convention adopted in the main part we use 
k = k 2 — K\ = 2k_, and the color factors are suppressed. The individual steps of the calculation 



are performed by a FORM-program [44 



(I) Firstly, the 7-algebra is evaluated. The result is represented by projections onto the matrices 

1, 7s, 7m, 7s7m, <V> 75<V • ( B - 8 ) 
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Here the only contribution is due to the axial vector part. The resulting expressions have been 
classified with the help of the pole-parts <//^'" Pn °f basic integrals, 



h=~\ 757, {2J& 

+ J32 ^(Pl + P2 pXa) ~ XpX<r{Pl + P2) ~ 9pXv%(Pl + P2) 



J Si 



31 



X tM (p lp Xp 2 +p 2 pXPl) ~ X p (j%Xp2+I%Xpi) 



(B.9) 



-2j£x% + J 20 x II x(p 1 +p 2 )}. 



In the calculation of the pole-parts of the individual integrals we apply the a-representation 
of the Feynman denominators and perform the momentum integrals. The exponentials are 
expanded into Taylor series as 



exp(iKkx) = 



(inkx) 7 



n=0 



ni 



Because of x 2 = only the first few terms of these expansions contribute. 
(II) The terms J^" pn in Eq. fljTg) are 

k x k»k v 



Jo 



'.33 



p.p. / d N k 

r 2 



q\q\k 2 

J D Q a^ a)K {(g Xl "P v {a) + g^ u P x {a) + g uX P^{ 



ITT 

+tn{S: x P fl (a)P u (a) + x> I P u (a)P x (a) + x" P x (a)P> 1 (a) 
D(a 



—ZK- 



[g^x" + g^x x + g uX x") 



,D(a) 



P x {a)x»x v + P^ajf^ + P v (a)x x 2P) + c^x^x"}., 



jpu 

'J □ ' 



'32 



qfqik 2 



2 

^ J D ae p{a);iK {(g^ + iK(P^(a)x v + P v {a)xf) +c 2 x fl x v }, 



J 31 



J 21 



p.p. / d N k^^-e ikiK = — z«x" / D n ae &p{ - a)K , 
qfqik 2 2e 

k,. 



p.p. / d N k^e ikiH 



ITT 

e 



D a5{l -ax- a 2 )e" p{a)K [P"{a) - -nx u D{a) 



J- 



20 



P.P. / dk 



ITT 

e 



Doa5(l — oc\ — «2)e 



ixP(a)n 



with the abbreviations 



D(a) = p\oli{1 - ai) +p\a 2 {l - a 2 ) - 2pip 2 aia 2 . 



(B.10) 



(b.ii; 



(B.12) 

(B.13) 
(B.14) 

(B.15) 
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The constants c x and c 2 in Eqs. (|B.11| , [B~T2|) remain unspecified here, since they do not contribute 
to the final result. Note that in the present case the two possible axial denominators 1/kx are 
canceled by corresponding factors in the numerator. For completeness we note that in integrals 
with axial denominators, like 



dk- 



1 



(p x — k) 2 (p 2 — k) 2 k 2 [kx 
we proceed as follows. The integrand is first rewritten as 



(B.16) 



kx 



kx 



kx 



(B.17) 



By applying the improved gauge prescription given by Leibbrandt and Mandelstam |45[ the 
integrals decompose into a part being calculated by Leibbrandt already, whereas the other part, 
important for the operator renormalization, can be calculated without difficulties. For diagrams 
containing more than one axial dominator similar decompositions hold, as e.g. 



— 1 — ikxK 



AkxK 



(kx) 2 N 
1 1 



Akxi 



(kx) [(k —p)x] N 



(kx) 2 N 
px N 



px N 



(k — p)x 

a ipxK 



1 .11 

+ JkxYH +m ~kx~H 

1 

kx 



(B.18) 



px N 



(k — p)x 

1 1 

kx 



+ e 



ipXK 



1 1 



(k — p)xpx N 



kx 



(B.19) 



Here N denotes a standard Feynman denominator. As a result for the pole-part of I x , Eq. ( |B . 71 ) 
we obtain 



r / \ / -N^V 

= -P.P.7 5 (7£)^j— 



D ae 



inxP(a) 



{5(1 



— a x — a 2 



2Px(a) - x(pi +P2) 



+ 



in(xpi)(xp 2 ) — 2xP(a) — in(xP(a)) 2 j 



(B.20) 



(III) Now some appropriate identities are applied to cast the expressions into a form containing 
the operators again. In our case these rules can be derived by partial integration. We use the 
following relations: 



D ae iKiP(a) [(iKxP(a)Y 
J D ae iK * p(a) iKx.P(a) 
J D ae iK£P{a) iKx Pl xp 2 



D ae iKiP{a) [5(l - a x - a 2 )(iKxP(a) - 3) + 6] (B.21) 
D ae iK * p{a) [5(1 - a x - a 2 ) - 2] (B.22) 



J D ae iKiP{a) 1 6(1 - a x - a 2 



xp 2 - —5 (a^ 

IK 



+—6(a 1 )5(a 2 ) (B.23) 

IK 



f D ae iK * p{a) 5(l - a x - a 2 )a 2 (xp x - xp 2 ) = [ D ae iK * p{a) — 6(1 - a x - a 2 )(l - 5(a x )) 

J J IK 



(B.24) 
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(IV) As the result of this calculation we obtain 



h 



-75(^7) I D ae i{K2 - Kl)Hpiai+P2a2) \5(a 1 )5(a 2 ) - 2} + finite terms 

e 4 J k 2 — K\ L > 



(B.25) 



Collecting all terms one gets 

1 rd N Pl d N p 2 



(rof{K 1 ,K 2 )s 



(2*0 



37V 



Xp 1 —tK 2 Xp2 



[-h(pi,P2, «1, «k) + A(Pl,P2, «2, «l)] 



+e i ™-™[-J 2 (p 2 ,p 1 ,K 2 ,K 1 ) +/ 2 (p 2 ,p 1)Kl ,K 2 )]}^(p 2 ) • 

(B.26) 



After inserting the explicit expressions and exploiting the symmetry relations of Eq. ( |B.25|) one 
obtains finally 



P.P.(TOf(K u K 2 )S 



d N pid N p 2 
(2tt) 8 



(B.27) 



x | e «(K2-Kl)i(piQl+P2a2) _|_ g-i(re2— Ki)2(piai+P2a2)| 

Here we have introduced the color factor Cp- Now we introduce new variables the 
possibility of changing a\ <-> a 2 and obtain 



^ dK' X dK' 2 I - ^ lC L P2 #G0l)75 (^7) V>(P2) 



X 



4(2tt) 4 7 1 z 7 (2tt) 8 

fgire'iSpi— in'2xp2 _|_ ^iK , 2xp 1 —iK' 1 xp2 



(B.28) 



-Dqo:- 



x 



«2 — «1 

- [«i + ai(«i - « 2 )])5(k' 2 - [«2 - Qt 2 («i - « 2 )]) 

- [«i - Qfi(ACi - K 2 )])8(k' 2 - [k 2 + Qf 2 («l - « 2 )]) 

Performing the momentum integration we finally end up with 



{5( ai )5(a 2 )-2} 



p.p ATO^kuk^S 



-J dK\dK' 2 -\^(K\x)j 5 (x^)i)(K' 2 x) +^(k' 2 x)j 5 (xj) ${k\x)) 

2ix 2 g 2 f 1 
x , — — / D a 



4(2tt) z 



K 2 — K\ 



C F {5( ai )5(a 2 ) -2} 



x 8{k'i - [k\ + ai(/ci - k 2 )])8{k' 2 - [k 2 - a 2 (^i ~ k 2 )\) 

+5(k\ - [m - a x {Ki - k 2 )])5(k' 2 - [k 2 + a 2 («i - k 2 )]) ■ (B.29) 

From this expression we can now read off the anomalous dimension, 

AK^( ai , a 2 ) = C F [5( ai )5(a 2 ) - 2] . (B.30) 

Finally we would like to add a remark concerning the Feynman rules. It is important to note 
that for non-local operators the involved field operators appear with unrestricted momentum 
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variables, e.g. A u (pi), so that we are not enforced to determine operator vertices explicitly. We 
can simply apply the standard Feynman rules if we have in mind that the non-local operators 
are in fact x-space operators so that for diagrams in momentum space there remain certain 
exponential functions as relic terms from the lacking ^-integration of the Fourier representation. 



C The [ ] + — prescription 

The anomalous dimensions naturally appear in a regularized form, Eqs. ( |6.57| - [6l60| ) and (|6.62j - 
S.65|) with the [ ] + -prescription given by ( |6.61| ). This prescription influences the subsequent 
relations, especially those of the evolution kernels V^(t,t' ,r). Let us consider the typical term 



ol\ + a 2 



5[t - t'(l -ai- a 2 ) + r(«i - a 2 )\C.l 



V sing (t, t', t) = / D a [6(atx) + 5{a 2 )\ 

which appears in V™t — AV^-t, KxP an d ^Kxt G - We introduce the variables 

«i = A£, a 2 = A(l - £) . 
The corresponding measures are 

/ D a [5(ai) + S(a 2 )} = C dX [ D£ = [ dX [ d£ [6(0 + 5(1 - 0] • 
J Jo J Jo Jo 



(C.2) 



(C.3) 



One obtains 



V sing (t,t',r) = J Di jT 1 dA[i] 5(t-t'(l-A) + rA(2e-l 

[5(T(X)-tX)+5(T(X)+tX)] 



dX 
"i dX 



[ [6(T(X) - rX) + 5(T(X) + rA) - 25(t - t')} 
Jo X 



(C.4) 



with T(A) = t — t'(l — A). The latter integral may now be rewritten in terms of the variable x 
and y, Eq. ( |7.36| ), yielding 



1 



2\t\ [y - x}_ 



[p{x,y) - p(x, 



2|r| 



p(x,y) p(x,y) 



[y - x], [y - x]_ 



(C.5) 



where the 



-prescription acts to the right. 



D Sample calculation of extended evolution kernels 

This Appendix is devoted to show how the extended evolution kernels are obtained. As an 
example we consider V e ^(t, t', r). We evaluate Eq. ( |7.34| ) 

Vext{t,t,i~) = dai da 2 K H (a 1 ,a 2 ) — 

JO JO J-oc Z71 

x (inxp + t)^ 1 exp {%Kxp + [t — (1 — a\ — a 2 )t' + r(a\ — a 2 )}} , (D.l) 
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with 



Because of the scaling relation 



-C F {S(a 1 )5(a 2 )+2}. 



r ext r'r' 



(D.2) 



(D.3) 



and reminding the general investigations in ||, we perform the calculations for r = 1 and 
< x,y < 1. The relation of the variables x and y to t,t' and r is given in Eqs. ( |7.36| ). The 
internal integral over nxp + yields 



The latter equation is distribution-valued in the sense of tempered distributions 
l.h.s of Eq. ( p.4|) is the Fourier-transform of the distribution 

2m z 



(DA) 



34 . The 



This distribution can be represented by 
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From Eqs. (|8TT5| a) and (|57T7| ) follows then Eq. ([TX4]). 

In terms of the variables x, y we solve 
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where 
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(D.9) 



J D aQ [x — (1 — a% — a 2 )y — a 2 ] 5(ai)5(a 2 ) 

Furthermore we observe that 

J D a<d(^ — [x — (1 — a\ — a 2 )y — 0:2]) = J D a<d(x — (1 — a,\ — a 2 )y — a 2 

Using these equations and taking into account the general structure of V ex 9 (x, y) we obtain finally 
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Eq. (p.lOj) also allows to derive the splitting function in the case of forward Compton scattering. 
It is obtained in the the limit r — > for t' > 0, with sign(y) = — sign(y) = 1. Here the 
variables x, y are expressed in terms of the variables (£,£') and r, Eqs. ( |7.36| ). In this limit 
x/y = x/y — > z = t/t' holds. Thus V e ^(t, t', r) is given by 



t'>0 



1 { ot 2 — Aft' 1 

— C F Q(z)Q(l -z)h + —tfs—j + 0(r) (d.ii; 



for small values of r, from which 
results. 
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